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Abstract 

Lagrangian reduction by stages is used to derive the Euler-Poincare equations for the nondis- 
sipative coupled motion and micromotion of complex fluids. We mainly treat perfect complex 
fluids (PCFs) whose order parameters are continuous material variables. These order param- 
eters may be regarded geometrically either as objects in a vector space, or as coset spaces of 
Lie symmetry groups with respect to subgroups that leave these objects invariant. Examples 
include liquid crystals, superfluids, Yang-Mills magnetofluids and spin-glasses. A Lie-Poisson 
Hamiltonian formulation of the dynamics for perfect complex fluids is obtained by Legendre 
transforming the Euler-Poincare formulation. These dynamics are also derived by using the 
Clebsch approach. In the Hamiltonian and Lagrangian formulations of perfect complex fluid 
dynamics Lie algebras containing two-cocycles arise as a characteristic feature. 

After discussing these geometrical formulations of the dynamics of perfect complex fluids, 
we give an example of how to introduce defects into the order parameter as imperfections (e.g., 
vortices) that carry their own momentum. The defects may move relative to the Lagrangian 
fluid material and thereby produce additional reactive forces and stresses. 
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1 Introduction 



Definitions. The hydrodynamic motion of a complex fluid depends on variables called order 
parameters that describe the macroscopic variations of the internal structure of the fluid parcels. 
These macroscopic variations may form observable patterns, as seen, for example, via the gradients 
of optical scattering properties in liquid crystals arising due to the spatially varying orientations 
of their molecules, as discussed in, e.g., Chandrasekhar [1992] and de Gennes and Prost [1993]. 
This micro- order of a complex fluid is described by an auxiliary macroscopic continuum field of 
geometrical objects associated with each fluid element and taking values in a vector space (or a 
manifold) called the order parameter space. The canonical example is the description of the local 
directional asymmetries of nematic liquid crystal molecules by a spatially and temporally varying 
macroscopic continuum field of unsigned unit vectors called "directors," see, e.g., Chandrasekhar 
[1992] and de Gennes and Prost [1993]. 

Thus, the presence of micro-order breaks the symmetry group O of the uniform fluid state 
to a subgroup V. This spontaneous symmetry breaking occurs in every phase transition of 
condensed matter into an ordered state. The symmetry subgroup V C O that remains is the 
isotropy group of whatever geometrical object appears (e.g., a vector, a spin, a director, etc.) when 
the symmetries of the uniform fluid state are broken to produce the micro-order. That is, the 
remaining symmetry subgroup V is the symmetry group of the micro-order. Equivalently, the 
order parameter may also be regarded as taking its values in the coset space C = O/V and its space 
and time variations may be represented by a space and time dependent curve in the Lie symmetry 
group O through the action of O on its coset space C. Thus, associating a geometrical object in 
a vector space (say, a director) with an order parameter is a way of visualizing the coset space C, 
see, e.g., Mermin [1979] and Volovick [1992] for physical examples. The symmetry group O of the 
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original uniform fluid state is called the order parameter group, or the broken symmetry. The 

order parameter dynamics that generates a curve parameterized by space and time in the broken 
symmetry group for a complex fluid is called its micromotion, although it refers to continuum 
properties at the coarse-grained macroscopic scale. 

Spatial and temporal variations in the micro-order are measured relative to a reference con- 
figuration. Let O act transitively from the right on a manifold M. Suppose the subgroup V C O 
leaves invariant an arbitrarily chosen reference point niQ € M, i.e., niQp = niQ , € V. The refer- 
ence point niQ then corresponds to the coset [e] = eV = "P of O/V, where e is the identity element 
of O. For another choice of reference point, say uiq = rriQh with h ^ O, the isotropy subgroup 
becomes the conjugate subgroup V ^ V = h~^Vh and the corresponding coset space O/V' is, 
thus, isomorphic to the original coset space O/V. Hence, the conjugacy equivalence classes of the 
order parameter coset spaces account for the arbitrariness in the choice of reference point niQ. One 
may think equally well of the order parameter group as acting either on a manifold M with "origin" 
mo, or on a coset space O/V, where V is the stabilizer of the reference point niQ. 

For example, let O be the group of proper orthogonal transformations 5*0 (3) acting transitively 
on directors in (unit vectors with ends identified) and choose mo to be the vertical director in an 
arbitrary reference frame; so mo is invariant under the 0(2) group of rotations around the vertical 
axis in that frame and reflections across its horizontal plane. This example applies to cylindrically 
symmetric nematic liquid crystal properties. The isotropy group V is 0(2), and the liquid crystal 
director may be represented in the coset space S'0(3)/0(2), which is isomorphic to S"^ /Z2, the unit 
sphere with diametrically opposite points identified, i.e., the real projective plane KP^ . Thus, 
after a reference configuration has been chosen, the micro-order of nematic liquid crystals may be 
represented equivalently as a space and time dependent curve in the broken symmetry group O 
acting on either the order parameter manifold M = RP^, or on the coset space S'0(3)/0(2) of the 
broken symmetry 5*0(3). 

In addition to its micromotion, the motion of a complex fluid also possesses the usual properties 
of classical fluid dynamics. In particular, the motion of a complex fluid involves the advection of 
thermodynamic state variables such as heat and mass, regarded as fluid properties taking values in 
a vector space V* . The complete dynamical equations for complex fluids must describe both their 
motion and their micromotion. In general, these two types of motion will be nonlinearly coupled 
to each other. This coupling typically arises because properties of both types of motion appear in 
the stress tensor governing the total momentum. 

This deflnition of complex fluid motion in terms of the continuum dynamics of both its order 
parameter and its usual fluid properties encompasses a wide range of models for complex fluid 
motion, including binary fluids, multi-phase fluids, polymeric materials, spin glasses, various other 
types of magnetic materials, superfluids and, of course, liquid crystals. The order parameter for 
each of these models provides a continuum (i.e., coarse-grained) description of the complex fluid's 
internal degrees of freedom, or micro-order. Here we will discuss mainly the case in which these 
order parameters are continuous material variables, that is, they are continuous functions 
carried along with the fluid parcels.^ Such media shall be called perfect complex fluids. This 
name (abbreviated PCF) is chosen to contrast with the perfect simple fluid (which has no internal 
micro-degrees of freedom) and to provide a geometrical basis for treating defects as imperfections in 

^In some cases, order parameters are determined from constraint relations that are Eulerian in nature, e.g., the 
volume fraction in two-phase interpenetrating flow as in Holm and Kupershmidt [1986] . This case will not be discussed 
here. 
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the order parameter field that can propagate, or move relative to the material labels. The relative 
motion of defects through the medium generally introduces additional reactive forces and sources 
of dissipation. We believe the dynamics of defects in complex fluids is best approached after first 
discussing the dynamics of perfect complex fluids, which is interesting in its own right. Here we 
shall concentrate on deriving the nonlinear dynamical equations for the ideal continuum motion 
and micromotion of perfect complex fluids. 

Once the equations for the nonlinear dynamics of their ideal (nondissipative) motion and micro- 
motion are established, dissipative processes must be included for most physical applications of 
complex fluids. By tradition, this is accomplished phenomenologically in these models, by intro- 
ducing kinetic coefficients, such as viscosity, mobility, thermal diffusivity, etc., so as to obey the 
requirements of the Clausius-Duhem relation that the entropy production rate be positive when the 
dynamics of all thermodynamic variables (including order parameters) are included, as in Dunn and 
Serrin [1985], Hohenberg and Halperin [1977]. Here we shall ignore dissipation entirely, trusting 
that it can be added later by using the standard phenomenological methods. 

For the case that the order parameter group O is the proper orthogonal group SO (3), a geo- 
metrical approach to complex fluids exists as part of the rational theory of continuum dynamics for 
materials with orientational internal degrees of freedom, such as liquid crystals. Rational theories 
of such complex media began with E. and F. Cosscrat [1909]. The Cosserat theories were recapit- 
ulated at various times by many different people. Sec, e.g., Kleman [1983] and Eringen [1997] for 
descriptions of recent developments and proposed applications of this approach for treating, e.g., 
liquid crystal dynamics in the tradition of the rational theory of continuum media. 

The present paper starts with the example of the Ericksen-Leslie theory of nematic liquid 
crystals and develops the geometrical framework for continuum theories of perfect complex fluids. 
In this geometrical framework, the motion and micromotion are nonlinearly and self-consistently 
coupled to one another by the composite actions of the diffeomorphisms and the order parameter 
group. The micromotion follows a curve in the order parameter group depending on time and 
material coordinate, and the motion is a time-dependent curve in the group of diffeomorphisms, 
which acts on the material coordinates of the fluid parcels to carry them from their reference 
configuration to their current positions. A feedback develops between the composite motion and 
micromotion, because the stress tensor affecting the velocity of the diffeomorphisms depends on 
the gradient of the order parameter. The mathematical basis for our development is the method 
of Lagrangian reduction by stages, due to Cendra, Marsden and Ratiu [1999]. As we shall see, 
obtaining the Euler-Poincare equations for perfect complex fluids requires two stages of Lagrangian 
reduction, first by the order parameter group and then by the diffeomorphism group. 

The main results in this paper. The Euler-Poincare approach provides a unified framework 
for modeling the dynamics of perfect complex fluids that preserves and extends the mathematical 
structure inherent in the dynamics of classical fluids and liquid crystals in the Eulerian description. 
This paper provides detailed derivations and discusses applications of the Eulcr-Lagrange equations, 
the Lagrange-Poincare equations, the Euler-Poincare equations, the Clebsch equations and the Lie- 
Poisson Hamiltonian structure of perfect complex fluids. The action principles for these equations 
appear at their various stages of transformation under reduction by symmetries. The theme of this 
paper is the reduction of degrees of freedom by transforming to variables in the action principles 
that are invariant under the symmetries of the Lagrangian. (This is the Lagrangian version of 
Poisson reduction on the Hamiltonian side.) The transformation is done in stages, not all at once. 
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for the additional perspective we hope it brings and to illustrate how Lagrangian reduction can 
be iterated in condensed matter applications. It may also be possible to impose nonholonomic 
constraints by stages, by transforming to variables that either respect the constraints, or appear in 
their specifications. 

The main new results for liquid crystals and PCFs in this paper are: 

1. Four action principles for PCF dynamics and their associated motion and micromotion equa- 
tions at the various stages of reduction - the Euler-Lagrange equations, the Lagrange-Poincare 
equations, the Euler-Poincare equations, and the Clebsch equations. 

2. The canonical and Lie-Poisson Hamiltonian formulations of these equations, their Lie alge- 
braic structures and the Poisson maps between them. 

3. The momentum conservation laws and Kelvin-Noether circulation theorems for these equa- 
tions. 

4. The reduced equations for one-dimensional dependence on either space or time, and the rela- 
tion of these reduced equations to the Euler- Poisson equations for the dynamics of generalized 
tops. 

5. A strategy for composing the dynamics of defects in a complex fluid with its underlying PCF 
dynamics. 

Outline. In Section ||, we develop these results for the motion and micromotion in the example 
of nematic liquid crystals, in forms that shall parallel the results of the general theory derived 
in the following section. The Euler-Lagrange equations from Hamilton's principle specialize to 
the Ericksen-Leslie equations, upon making the appropriate choices of the kinetic and potential 
energies. The Lagrange-Poincare equations and Euler-Poincare equations that follow from applying 
two successive stages of Lagrangian reduction of Hamilton's principle with respect to its symmetries 
provide geometrical variants and generalizations of the Ericksen-Leslie equations. 

In Section ^, we perform two successive stages of Lagrangian symmetry reduction to derive 
first the Lagrange-Poincare equations and then the Euler-Poincare equations for PCFs with an 
arbitrary order parameter group. These Euler-Poincare equations are Legendre-transformed to 
their Lie-Poisson Hamiltonian form. Their derivation by the Clebsch approach is also given. 

Finally, in Section ^ we discuss a strategy of including defect dynamics into the theory of PCF 
dynamics. This strategy introduces an additional set of fluid variables that describe the motion 
of the defects relative to the material coordinates of the PCF. We shall illustrate this strategy in 
the example of rotating superfluid "^He, in which U{1) is the broken symmetry and the defects are 
quantum vortices. 

2 The example of liquid crystals 
2.1 Background for liquid crystals 

We begin with a hands-on example that illustrates the utility of the ideas we shall develop in this 
paper. Liquid crystals provide a ubiquitous application that embodies these ideas and supplies a 
guide for developing them. 
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Liquid crystals are the prototype for complex fluids. For extensive reviews, see Chandrasekhar 
[1992] and de Gennes and Prost [1993]. An orientational order parameter for a molecule of arbitrary 
shape is given in Chandrasekhar [1992], p. 40, as 

S = ^(^3x^X-Id®Ict^, (2.1) 

where ( • ) is a statistical average and x ^ 5*0 (3) is a rotation that specifies the local molecular 
orientation relative to a fixed reference frame. Thus, in index notation, 

SklKL = -^(jiXkK XiL - SmSklJ ■ (2.2) 

This order parameter is traceless in both pairs of its indices, since = X^^- 

For cylindrically shaped molecules, (e.g., nematics, cholesterics or smectics^ we may choose 
the 3-axis, say, as the reference axis of symmetry, i.e., choose K = 3 = L and set 

nk = Xk3, so that \nf = xskXks = 1 , (2.3) 

and S becomes, 

Ski33 = Ski = ^(3nkni-6ki^- (2.4) 

Note: The order parameter S does not distinguish between n and — n. Physically, S may be 
regarded as the quadrupole moment of a local molecular charge distribution. For a clear description 
of the use of this order parameter in assessing nematic order-disorder phase transitions using the 
modern theory of critical phenomena, see Lammert et al. [1995]. Instead of considering such phase 
transitions, here we shall be interested in the continuum dynamics associated with the interaction 
of this order parameter with material deformations. 

In passing to a continuum mechanics description, one replaces statistical averages by a local 
space and time dependent unit vector, or "director" n(x, i). Then, the continuum order parameter 
S corresponding to the statistical quantity S in (|2.4| ) is the symmetric traceless tensor. 

Ski (n) = ^(3nkni- 6ki^ , (2.5) 

which satisfies S{n) = S{—n) and admits n as an eigenvector, 5 • n = n. The hydrodynamic tensor 
order parameter S represents the deviation from isotropy of any convenient tensor property of the 
medium. For example, the residual dielectric and diamagnetic energy densities of a nematic liquid 
crystal due to anisotropy may be expressed in terms of the tensor order parameter S as, 

^ E • 5(n) • E and ^ B • S{n) ■ B , (2.6) 

for (external) electric and magnetic fields E and B, respectively. (For simplicity, we neglect any de- 
pendence of the electric and magnetic polarizabilities of the medium on the gradients Vn, although 
this possibility is allowed.) 



^Smectics form layers, so besides their director orientation they have an additional order parameter for their 
broken translational symmetry. 
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2.2 Four action principles for perfect liquid crystals 

Liquid crystal action. The standard equations for the continuum dynamics of liquid crystals 
without defects are the Ericksen-Leslie equations, due to Ericksen [1960, 1961], Leslie [1966, 1968], 
and reviewed, e.g., in Leslie [1979]. These equations express the dynamics of the director n and 
may be derived from an action principle (^5 = with action <S = j dtL in the following class. 



S = J dtj d^X >C(x, J, n, h, Vn) , 



(2.7) 



with notation x = (9x(X, t)/9t, so that overdot denotes material time derivative, J = det((?x/9X) 
and Vn(X, t) has spatial components given by the chain rule expression. 



V,;n 



\dxJ 



-1 



dn OXa d 



AiA 



OXa dxi BXa 



n(X, t) = n i 



(2.8) 



Note that the coupling between the fluid dynamics x(X, t) and director dynamics n(X, t) occurs in 
the Lagrangian density C through the inverse deformation gradient, (5x/9X)~^, via the chain 
rule expression above for Vin(X,t). 

Varying the action S in the fields x and n at fixed material position X and time t gives 



5S 



d^Xl 5x 



dC 

dxjj 



d dC 8 
+ Jt^— - J- 



dxp dJ 



dXr. 



+ 6n ■ 



. \ dn ) dn 



dXr, 



J ^ n,p 

onm 



dC 
dnrr, 



(2.9) 



with natural (homogeneous) conditions, expressing continuity of 

dC , . dC 
— and n„i- , 



(2.10) 



on a boundary, or at a material interface, whose unit normal has spatial Cartesian components n^. 
(We always sum over repeated indices.) 

Action principle #1 — Euler-Lagrange equations. Stationarity of the action, SS = 0, i.e., 
Hamilton's principle, thus yields the following Euler-Lagrange equations in the Lagrangian 
fluid description. 



Sxp : 



dC\ B dC a ( . dC 

+ Jt^— - J- — I J n 



dx. 



dxp dJ dx . 



dnm 



0, 



<5n- (—\- — +J-^(j-^^^\-Q 
\dn) dn dxm \ dn^m J 



(2.11) 
(2.12) 



The material volume element J = det{dx/dX.) satisfies an auxiliary kinematic equation, obtained 
from its definition. 



J = J 



dXA dxi 
dxi dXA 



(2.13) 
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Imposing constant J, say J = 1, gives incompressible flow. The last term in equation is the 

divergence of the Ericksen stress tensor arising due to the dependence of the potential energy of 
the medium on the strain Vn, see Ericksen [1960, 1961]. 
With the proper choice of Lagrangian density, namely, 

£ = ||x|2 + ^|h|2 +p{J - 1) + g(|n|2 - 1) - F(n, Vn) , (2.14) 



where p and / are material constants one finds that the Euler-Lagrange equations ( p. 11 ) and (2.12) 



produce the Ericksen- Leslie equations for incompressible liquid crystal flow, Chandrasekhar 
[1992] and de Gennes and Prost [1993], 

dF d dF 

/n-2gn + h = 0, with h = — - — . (2.16) 

an oxj on J 



In the Lagrangian density ( p. 14 ), the Lagrange multipliers p and q enforce the incompressibility 



condition J = 1 and the director normalization condition |n|^ = 1, respectively. The standard 
choice for the function -F(n, Vn) is the Oseen-Zocher-Frank Helmholtz free energy density, as 
discussed, e.g., in Chandrasekhar [1992] and de Gennes and Prost [1993], namely,^ 

F(n,Vn) = A:2(n • V X n) + — (V • n)2 
chirality ^^^^^ 

+ :^(n.Vxn)2 + M|nx Vxn|2 
2 ^ ' 2 ' ' 

twist bend 
+ ^E.5(n).E + ^B.5(n).B, (2.17) 

with /c2 = for nematics (but nonzero for cholesterics) and n • V x n = for smectics. Since 
smectics form layers that break translational symmetry, their order parameter group may be taken 
as the Euclidean group E{3), or perhaps as 50(3) x U{1) in simple cases. We shall concentrate our 
attention on nematics in this section. However, the general theory for arbitrary order parameter 
groups developed in the next section would also encompass smectics. 

Key concepts: Material angular frequency and spatial strain of rotation. The rate of 
rotation of a director n in the rest frame of the fluid material element that carries it is given by 

i/ = nxh. (2.18) 

In terms of this material angular frequency v (which is orthogonal to the director n), the 
dynamical equation (|2.16| ) becomes simply, 

/j> = hxn. (2.19) 



3 



This notation treats spatial and director compone nts on the same (Cartesian) basis. Later, we will distinguish 



between these types of components, see, e.g., equation (2.82) 
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Likewise, the amount by which a specified director field n(x) rotates under an infinitesimal spatial 
displacement from Xi to xi + dxi is given by, 

7j = n X Vj n or, j ■ dx = -Jidxi = nx dn. (2.20) 

Remarks on kinematics. An orthogonal matrix with components XkL is associated with an 
orthonormal frame of unit vectors (1, m, n) whose components {Ik^f^ki^k) = {XkiiXk2iXkz) satisfy 

^jkl{dxX~^)lk = (1 X dl + m X dm + n X dn)j , (2.21) 

which is 0(3) right invariant. Accordingly, the spatial rotational strain, with components 

(n X dn)j = ejki{dxi2,Xkz) , (2.22) 

and the material angular velocity, 

(n xh)j = ejkiixBXks) , (2.23) 

are each elements of TS'^ /Z2, i.e., they are tangent to the sphere |np = 1 and 0(2) right invariant 
under 5*0(2) rotations about n and under reflections n — > — n. Also, n x dn and n x n are 0{2) right 
invariant elements of the Lie algebra, so(3) ~ TSO{3)/SO{3), so they are elements of so(3)/0(2). 
Being orthogonal to n, their evolution preserves |np. Therefore, |np = 1 may be taken as an initial 
condition, and n may be reconstructed from h = v x n. 

The spatial rotational strain n x dn = 7 • dx satisfies the kinematic relation, 

(7 ■ dx) • = (2i^ X 7i + Vj u) dxi , (2.24) 

obtained from its definition. This kinematic relation involves only 7^, v and Xj, which suggests 
transforming variables from (n, n, Vn) to 7^ and u. (Note that 7j and u do not distinguish between 
n and — n.) Each component of 7^ is orthogonal to the director n. That is, n.7j = 0, for i = 1,2, 3. 
Using this fact and |np = 1 gives the relations, 

Vi n = - n X 7j and V x n = - n tr7 + n"* 7^ , (2.25) 

or, equivalently, in index notation, 

V^n" = -e"'''=nSf and e^^'^ V^n^^ = -n^m + 7^ > (2-26) 

where lower indices are spatial and upper indices denote director components. Consequently, we 
have the following useful transformation formulas. 



Vj n • Vj n = 7j • -fj , Vj n x Vj n = 7^ x 7^ and n • V x n = — tr7 . 



(2.27) 
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Connection, curvature, singularities and topological indices. The Eulerian curl of the 
spatial rotational strain 7^ yields the remarkable relation, 

= - + 27^ X 7^. = , (2.28) 

where Bij vanishes, provided the vector field n is continuous. In geometrical terms, vanishing of Bij 
is the Maurer-Cartan relation for the flat connection one-form (or, left-invariant Cartan one-form) 
fidxi = n X dn, as discussed in, e.g., Flanders [1989]. That is, the curvature two-form Bijdxi A dxj 
vanishes in the absence of singularities (disclinations) in the director field n(x). Thus, By may 
be regarded as the areal density of disclinations in a nematic liquid crystal. In Section we shall 
discuss how to proceed when defects exist and the disclination density Bij does not vanish. 

A second remarkable formula also stems from the curl— 7 equation ( p. 28 ), which implies the 
relation, 

^dxi Adxi = e'^^'^^^dxi Adxi = e'^'^^dn^ A dn' , (2.29) 

OXi OXi OXj 

for a continuous director field. Contracting with and taking the exterior derivative of this 
relation implies 

d f d-f\ 



dxi 



fn°-^ jdxi Adxj A dxk = e^^^^dn" A dn^A dn" = det(Vn)d^x = detij)d^x . (2.30) 



By |np = 1, these expressions must vanish, when the director field has no singularities. Thus, in 
the absence of singularities, there is a vector v, for which 

ekijn"-^ = ekije^^'n" = euijVj^i = (curl v)fc , (2.31) 

or, equivalently, 

e'^'^^n'' Vn^ x Vn^ = curl v . (2.32) 

Up to an inessential multiplicative factor, this formula is the well known Mermin-Ho relation 
between three-dimensional superfluid texture n and vorticity curl v in dipole-locked ^He — A, found 
in Mermin and Ho [1976] . Since n is a unit vector, a simple calculation in spherical polar coordinates 
transforms this relation to 

Vn^ X Vn^ = V(/> X Vcos6' = curl v , (2.33) 

for n = (sin^cosc/), sm.9sm.(f), cos9)^ in terms of the polar angle 9 and azimuthal angle (f> on the 
sphere. Hence, 

curl v • dS = e^'^'n'^dn^ A dn" = # A dcos9 , (2.34) 

which is the area element on the sphere. Thus, by Stokes' theorem for a given closed curve C in 
with unit tangent vector n(x), the integral 

27rWr = (f vdx= If e^^'n^-^l dxk A 7f dxi , (2.35) 
Jc JJs 
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is the (signed) area swept out by n on a portion of the unit sphere upon traversing the curve C 
which is the boundary of the surface S in three dimensions. This is also equal to the writhe of the 
curve C, see Fuller [1978], as cited in Kamien [1998]. Thus, if the field lines of the three dimensional 
vector field n(x) form a closed curve C, the writhe of this curve is given by the area integral in 



equation ( |2.35| ) over a surface whose boundary is C. Equation ( ^.35| ) is the second remarkable 
formula implied by the curl— 7 equation ( 2.28| ). 



Finally, the integral of v • curl v over three dimensional space gives a number, 

N = ^ vcurlvd^x, (2.36) 



4tt 



53 



which takes integer values and is called the Hopf index or degree of mapping for the map n(x) : 
^ S'^. The Hopf degree of mapping A'^ counts the number of times the map n(x) : ^ S'^ covers 
the unit sphere, as discussed in, e.g., Flanders [1989]. Equivalently, the Hopf degree of mapping 
N counts the number of linkages of the three dimensional divergenceless vector field curlv with 
itself. See Mineev [1980] for the physical interpretation of this formula in terms of disclinations 
in nematic liquid crystal experiments. See also Trebin [1982], Kleman [1983], Kleman [1989] and 
references therein for additional discussions of the differential geometry of defects in liquid crystal 
physics. 

The quantity e^^'^n^Vn'' x Vn'^ in equation ( 2.33| ) appears in many differential geometric contexts 



in physics: in the Mermin-Ho relation for the vorticity of superfiuid ^He—A in terms of the "texture" 
n due to Mermin and Ho [1976]; in the Skyrmion Lagrangian, as discussed in, e.g., Trebin [1982]; 
as the n— field topological Wess-Zumino term in the 0(3) nonlinear sigma model, discussed in, 
e.g., Yabu and Kuratsuji [1999], Tsurumaru [1999]; and as the instanton number density in coset 
models. For references to the last topic see Coquereaux and Jadcyk [1994] and for a discussion 
of the associated Poisson-Lie models, see Stern [1999]. This same term also produces forces on 
vortices in ^He — A{B) as found in Hall [1985] and in ferromagnets, discussed in Kuratsuji and 
Yabu [1998]. Also, we have seen that this term also allows the calculation of the "writhe" of a 
self-avoiding closed loop of, say, DNA, or some other polymer filament with unit tangent vector 
n, see, e.g.. Fuller [1978], Goldstein et al. [1998], Kamien [1998]. See also Klapper [1996] and 
Goriely and Tabor [1997], for additional recent discussions of the differential geometric writhe of 
DNA conformations. Finally, we have identified the relation of the Hopf degree of mapping to the 
linkage number of the writhe flux, curlv. Perhaps surprisingly, in the case of liquid crystals, it 
will turn out that this topological index, the writhe of a fiuid loop, and the flux of curvature Bij 
through a fiuid surface (which measures the presence of singularities in the director field) are are 
all three created nonabelian PCF dynamics. 



Action principle ^2 — Lagrange-Poincare equations. Intrigued by the geometrical proper- 
ties of the variables 7j and v, we shall seek to express the Ericksen-Leslie equations solely in terms of 
these variables, by applying symmetry reduction methods to Hamilton's principle for liquid crystal 
dynamics. The trace of the spatial rotational strain 7 gives 

7^ = e'^^Vv^n" tr7 = 7;;j = -n-Vxn, (2.37) 



which figures in the energies of chirality and twist in equation ( p. 17 ) for F(n, Vn). 
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We shall now drop the anisotropic dielectric and diamagnetic terms, which exert torques on the 
director angular momentum due to external E and B fields given by, 

hEM X n = ■ n)E X n + ^^^iB -11)6x11, (2.38) 

but which do not contribute in the motion equation. Upon neglecting these torques, which may 
always be restored later (see the Appendix), the remaining Oseen- Frank free energy F{n, Vn) may 
be expressed (modulo a divergence term that will not contribute to the Euler-Lagrange equations) 
entirely as a function of a subset of the invariants of 7, namely, tr7, tr(7 + 7"^)^ and tr(7 — 7^)^^, 
with the help of the identities below, obtained using |np = 1. In the notation of Eringen [1997], 
one has 

h = (V X nf = [n • (V X n)]^ + [n x (V x n)]^ , (2.39) 

I2 = n^jn^j = (V - n)2 + V X n) V (^n^n^j - n'n^j^ (2.40) 

\ 2 

V X n) + n-'-n* 



h = [(n • V)n)]2 = [n X (V X n)]^ . (2.41) 
Geometrically, the divergence term 

(n^72^j - n'n^j^ = 2kiK2 , (2.42) 

appearing in equation ( p. 40 ) for I2 is (twice) the Gaussian curvature of the surface S whose normal 
is n(x). This equation follows from Rodrigues' formula 

(in = -Kdx, (2.43) 

for the principle directions of the surface S. Such a surface exists globally, provided n - curln = 0, 
see Weatherburn [1974], as cited in Kleman [1983]. Straightforward calculations give the following 
linear relations, see Eringen [1997] 

tr(7|) = -^/i + /2+ ^/3, tr(7i) = -i(/i+/3), (2.44) 
where 75 and 7^ are the symmetric and antisymmetric parts of 7, 

75 = ^(7 + 7^), and 7^ = 1(7-7^). (2.45) 

Hence, modulo the divergence term in ( |2.40| ), the Oseen- Frank free energy in ( 2.17| ) may be written 
in terms of the invariants of 7 as F{ tr7, tr(7|), tr(7^) ). 

The reduced Lagrange-Poincare action. By transforming to the rotational strain 7^ and 
frequency v, the action S for liquid crystals in ( p.7| ) may thus be reduced, that is, rewritten in 
fewer variables, by defining, 

S = jMjSxa^.J,u.,J, (2.46) 
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where the reduced variables u and 7^ are perpendicular to the director n. 

A calculation using their definitions shows that the variations of J and v satisfy 

OXa d6xi . \. „ / . N 

oJ = J— — , Oi/ = X c)nj — 2i/ X X on] , 

(JXi OJ\A 



(2.47) 



in terms of the variational quantity n x 6n. Likewise, we calculate the variation of the rotational 
strain 7^ as 



6j, 



J 



J 



n x on) 



OXa L dx 
27m X (n X 5n) - 7^ J 



^1 d 



QXa 



J- dXj: 

OXm. 



(2.48) 



This variational expression for 7^ may be rearranged into the more suggestive form, 

^Im + Iv^m fep = 2(n X (^n) X 7„ + Vm(n x (5n) , (2.49) 
in which its similarity with the kinematic relation for 7 in ( 2.24| ) is more readily apparent, especially 



when the kinematic equations (2.13) and (2.24) are rewritten in the notation of equation ( |2.8| ) as 

j = JVi Xi and 7^ + 7pVm Xp = 2i/ x 7^ + Vm . (2.50) 

The variation of the Lagrange- Poincare action S in ( p. 461 ) for liquid crystals in the Lagrangian fluid 
description may now be rewritten as 



bS = - 
+ (n X 5n 



dt / d-^Xi 6x 



dC \ J d dC J d 

dxp J dxp dJ dx„ 



J 7p 



din 



dc ^ d 

OU OXm 



J 



_i dC 

dim 



'^Im X 



dim 



+ dt d-^Xl 5x. 



+ 



d 



BXa 



dC 

'dxp 



dC \ 



+ (nx5„).- 

dC 



dXA . , ^ . 

dxm ^"^H-^aj"^"^ 97 



dXA 

dXrri 



n X 6n) 



dC 

dim 



(2.51) 



Consequently, the action principle 55 = yields the following Lagrange- Poincare equations for 
liquid crystals. 



5xp : 



dC\ d dC d 
+ Jt^— - J- 



dx. 



dxp dJ 



dXr. 



j^i dC 

J Iz 



dC 



n X (5n : — — — 2i/ x — h J 



d 



du J 



dv 



dXr. 



dlr, 

_i dC 



with natural boundary conditions, cf. equation ( 2.10| ), 



dC/dJ = and 



J 



dC 
dim 



dir. 



0, 



0, 



27m X 



dC 

dim 



0, 



(2.52) 



(2.53) 
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on a boundary, or material interface, whose unit normal has spatial Cartesian components hm- 
These are the dynamical boundary conditions for liquid crystal motion. These conditions 
ensure that the fluid pressure dC/dJ and the normal stress are continuous across a fluid interface. 
The term 7^ • dC/d^^y^ contributes to the stress tensor of the complex fluid and arises from the 
dependence of its free energy upon the rotational strain 7j. Again J = 1 for incompressible flow, as 
imposed by the Lagrange multiplier p, the fluid pressure. (The Lagrange multiplier q is no longer 
necessary.) 

Upon specializing to the Oseen- Frank-type Lagrangian (|2.14| ), stationarity of the action 5(x, J, 7) 
in equation ( |2.46| ) recovers the Ericksen-Leslie equations (|2.15| ) written in these variables. This is 
an example of Lagrangian reduction - from the variables (n, n, Vn) to {v, 7). The equations in 
fewer variables that result from the first stage of Lagrangian reduction were named the Lagrange- 
Poincare equations in Marsden and Scheurle [1995] , Cendra, Marsden and Ratiu [1999] . One of the 
goals of this paper is to characterize the geometrical properties of such equations for PCFs, using 
the mathematical framework of Lagrangian reduction by stages established in Cendra, Marsden and 
Ratiu [1999]. From this viewpoint, the first stage of the Lagrangian reduction for liquid crystals is 
finished. 



Action principle — Euler-Poincare equations. The Lagrange-Poincare equations (2.52) 



and their kinematic relations ( ^.50 ) are still expressed in the Lagrangian fluid description. We shall 



pass to the Eulerian fluid description by applying a second stage of Lagrangian reduction, now 
defined by the following right actions of the diffeomorphism group, 

u{x,t) = x{X,t)g-\t), 
D{x,t)d^x = d^Xg-^{t), 

u{x,t) = u{X,t)g-\t), (2.54) 

7,(x,t)dx, = (7.(X,t)^^gl^dX^)5-^(t), 

where the right action denoted x{X,t) = Xg{t) defines the fluid motion as following a time- 
dependent curve g{t) in the diffeomorphism group as it acts on the reference configuration of the 
fluid with coordinate X. In the Eulerian fluid description, the action S in Hamilton's principle 
6S = IS written as 



dt J d^x £{u,D,u,-f), (2.55) 

in terms of the Lagrangian density i given by 

e{u, D, u, f)d^x = C {±g-\t), Jg-^ (t), ug-\t) , -fg-^ (t)) (^d^Xg-^ (t)) . (2.56) 
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The variations of these Eulerian fluid quantities are computed from their definitions to be, 



drjj drjj 



dUj = + Uk 

5D 
5v 



dt dxk 
dDrjj 



dxh ' 



dxj 



as as 

dt dxm 



2u X "E — rj„ 



dv 

dXrn. 



27m X ^ - Vk 



dTn 

dxi 



Ik 



drjk 

dXrn. 



(2.57) 
(2.58) 
(2.59) 
(2.60) 



where S(x, t) = (n X 5n)(X, t)(7 ^{i) andrj = 5gg ^{t). One may compare these Eulerian variations 



with the Lagrangian variations in (2.47) and ( |2.48| ). See also Holm, Marsden and Ratiu [1998] for 
more discussion of such constrained variations in Eulerian fluid dynamics. Here the quantities D 
and 7^ satisfy the following Eulerian kinematic equations, also obtained from their deflnitions, 
cf. the Lagrangian kinematic relations ( ^.50 ), 



dD_ 

~dt 

dim 
dt 



dD 



Ui 



dxj 



du 

dXjy, 



+ 2 1/ X 7,^ 



dim 

dxk 



Ik 



duk 

dXrn. 



(2.61) 
(2.62) 



Note the similarity in form between the Eulerian variations of D and 7^, and their corresponding 
kinematic equations. This similarity arises because both the variations and the evolution equations 
for these quantities are obtained as inflnitesimal Lie group actions. 



Remark. The kinematic formula ( ^.62 ) for the evolution of 7^ immediately implies the following 
'y— circulation theorem for the spatial rotational strain. 



dt 



Im d-^r. 



c{u) 



2u X 7^ dXr, 



c{u) 



(2.63) 



Thus, the circulation of 7 around a loop c(n) moving with the fluid is only conserved when 7^ x v 
is a gradient. Otherwise, the curl of 7^ x v generates circulation of 7 around fluid loops. By 



Stokes' theorem and equation (|2.29D we have 

Im^-^m — / 
c(u) Js(u) 



dXj A dXr^ 

OXi 



Js{u) 



e'^^'dn^ A dn'' 



(2.64) 



for a surface S{u) whose boundary is the fluid loop c{u). Consequently, the 7— circulation theorem 
( 2.63| ) implies that a nonvanishing curl of 7^ x u generates a time-changing flux of e^^'^Vn'' x Vn^ 
for a = 1,2, 3, through those surfaces whose boundaries move with the fluid. 



The nonabelian 2-cocycle terms create withe and linkages. We use the kinematic equation 
(|2.62| ) for gamma and the deflnition of u to calculate the evolution of curlv ■ dS = n"e"^'^dn'' A dn^ 
as 



(^dt + £u) (curlv • fiS) = 2di>^ A dn'' = (i(n • 2v x 7^ dxm) = d{iY • a.diTfm dxm) 
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This means the nonabehan aspect of the generahzed 2-cocycle (the \7v in the gamma equation) 
creates writhe in the boundary of a surface moving with the fluid under the gamma-evolution. 
That is, by Stokes law, 



— (p V • ax 

(^t Jc{u) 



n • 2i/ X 7^ dxm 



n • adj/Yj^ dxm , 



' c(u) J c{u) 

and the writhe of a fluid loop is not preserved. Likewise, we calculate 

(5t + £u) (v • dx A d{v ■ dx)) / exact form . 

Therefore, we find 

d 



dt 



V • curl V d X / . 



Thus, perhaps surprisingly, the linkage number, the Hopf degree on mapping in equation (2.36) 
is also not preserved by the gamma-evolution. So the nonabelian generalized 2-cocycle term in 
the gamma equation also creates linkages in the Mermin-Ho quantity, curlv. Conversely, the 
preservation of these quantities is an Abelian characteristic. 



Euler-Poincare action variations. We compute the variation of the action ( p.55[) in Eulerian 
variables at fixed time t and spatial position x as, 



6S 



h 


fd^x 


h 


fd^xf^ 



-^6uj + —6D + — ■6u + 



+ I] 



d_M_ 
'diJu 



+ 



d 



dXr 



6D 

d_6e_ 
dt 6uj 

6£ dv 
5u dxj 



5v 

61 duk 



d 



6uk dxj dxk \5uj 
dim ^ d 



Uk) + D 



a 



dxj 6D 



Sir. 



dxj dxr. 



li 



Sin 



d f5i 51 

I Um H 



dXm ^SU 



+ 



d_ 
di 



Vj 



Sim 

51 



+ 2i/ X — + 27„ X — 
5v ™ (57^ 



5u, 



5D 



li 



51 

Sim 



5£_ 
6^ 



^ , J 51 
+ 2J • ( -— lim + - — 
5u ^7„ 



(2.65) 



where we have substituted the variational expressions ( 2.57 )-( 2^ ) and integrated by parts. 

The dynamical equations resulting from Hamilton's principle 5S = Q are obtained by requiring 
the coefficients of the arbitrary variations r]j and S to vanish. We assume these variations them- 
selves vanish at the temporal endpoints and we defer discussing the boundary terms for a moment. 
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Hence, we obtain the Euler-Poincare equations for liquid crystals, 

7? - -— = -—^-—{—uA + D—— (2 66) 



dt 6uj 6uk dxj dxk V 5uj ) dxj 5D 



+ 7r-|7,- 



du dXj 9xj dXm \ ^1 



d bi d , 

^ : — — = - 1 —Um + 



dtdu dxm^Sv 5 J, 

+ 21/ X — + 27^ X -— . (2.67) 

Equation (^^6^ ) agrees with the Lagrangian version of the micromotion equation governing u in 
(|2.52| ). To see this, it is helpful to recall that bljbv is an Eulerian density, so a Jacobian is involved 
in the transformation to the Lagrangian version. 

Momentum conservation. In momentum conservation form, the liquid crystal motion equation 
( 2.66| ) in the Eulerian fluid description becomes 



for simple algebraic dependence of the Lagrangian i on u, D, i> and 7^. This momentum conser- 
vation law is in agreement with the direct passage to Eulerian coordinates of the motion equation 
( |2.52| ) in the Lagrangian fluid description. For this transformation, it is helpful to recognize from 
equation ( p.56| ) that I = {CJ^^)g{t)~^ implies, by the chain rule, 

^g{t)-^ = I - D— , Since D{x,t) = J-\X,t)g{t)-K (2.69) 



The momentum conservation law ( 2.68| ) acquires additional terms, if the Lagrangian £ also depends 
on gradients of u, D , v and 7^. 

Noether's theorem. Noether's theorem associates conservation laws to continuous symmetries 
of Hamilton's principle. See, e.g., Olver [1993] for a clear discussion of the classical theory and 
Jackiw and Manton [1980] for its applications in gauge theories. The momentum conservation 
equation ( p.68| ) also emerges from Noether's theorem, since the action S in equation ( p.55| ) admits 
spatial translations, that is, since this action is invariant under the transformations, 

Xj — > x'j = Xj + r]j(x, t) with rjj = Cj , (2.70) 

for constants Cj, with j = 1,2,3. To see how equation ( p.68|) emerges from Noether's theorem, 
simply add the term d{ir]j)/dxj (arising from transformations of the spatial coordinate) to the 
endpoint and boundary terms in the variational formula ( p. 65 ) arising from variations at fixed time 
t and spatial position x, then specialize to r]j = Cj. 
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Kelvin-Noether circulation theorem for liquid crystals. Rearranging the motion equation 
dPI) and using the continuity equation for D in (|2.61| ) gives the Kelvin-Noether circulation 
theorem, cf. Holm, Marsden and Ratiu [1998], 



d_ r 1 6i 
dt Jc{u) D5uj ^ 



c(u) 



D 



51 , 51 
du+ — 
5v (57^ 



d 



dXr. 



(2.71) 



Hence, stresses in the director field of a liquid crystal and gradients in its angular frequency v 
and rotational strain 7 can generate fluid circulation. Equivalently, by Stokes' theorem, these 
gradients of liquid crystal properties can generate vorticity, defined as u = curl (Z)~^(5£/(5u). For 
incompressible flows, we may set D = 1 va. these equations and write the vorticity dynamics as. 



dt 



+ (u • V)uji - {u> ■ V)ui 



Vv"" X V 



51 



~\~ (^ijk 



5v"- 
_d d_ 



+ V7r„ X V 



5£ 



Ik • 



(2.72) 



Thus, spatial gradients in the director angular frequency v> and rotational strain 7 are sources of 
the fluid vorticity in a liquid crystal. 



2.3 Hamiltonian dynamics of perfect liquid crystals 

The Euler-Lagrange-Poincare formulation of liquid crystal dynamics obtained so far allows passage 
to the corresponding Hamiltonian formulation via the following Legendre transformation of the 
reduced Lagrangian i in the velocities u and u, in the Eulerian fluid description, 



nii 



5^ 
5ui 



h{m, D, (J, 7„) = rriiUi + a • v — £(u, D, v, 7^ 



Accordingly, one computes the derivatives of h as 

5h 5h 5h 

5^ 



5h 



5£ 



Ui 



V . 



(2.73) 



(2.74) 



5rai 5cT 5D 5D ^7^ ^Im 

Consequently, the Euler-Poincare equations ( p.66| ) - (2.67) and the auxiliary kinematic equations 



(2.61) - (2.62) for liquid crystal dynamics in the Eulerian description imply the following equations, 
for the Legendre-transformed variables, (m, Z), cr, 7), 



drrii 



dD 
'dt 
da_ 
'dt 

dt 



d 5h 



d 



dxi 5m j dxj 



+ 



djj_ 

dxi 



d 

dxj 



dxj 



D- 



5h 
■ (57, 
5h 



_d_ 

dxi 



5h 

'' 5mj 



dxj 5D 



5h 
<57,- 



d 5h 
dxj 6 a 



5mj 
5h 



5h_ 

5mj 5^j 



5h 

2(7 X — 

da 



27 ,- X — 



5h 
5-fj 



1 



d 5h 

dxi 6m j 



dli\ 5h 



dxj J 5m j 



+ 



d 6h 
dxi 5a 



5h_ 



(2.75) 
(2.76) 
(2.77) 
(2.78) 
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These equations are Hamiltonian. That is, they may be expressed in the form 

- = {z,M = b.-, 

where z € (m, D^a^'y) and the Hamiltonian matrix b defines the Poisson bracket 

5h 



(2.79) 



{/>} 



Oz dz 



(2.80) 



which is bihnear, skew symmetric and satisfies the Jacobi identity, 

{f,{g,h]]+c.p.{f,g,h)=Q. 
Assembhng the hquid crystal equations ( p. 75]) - ( 2. 781) into the Hamiltonian form ( 2.79 ) gives, 





mi 




mjdi + djmi 


Dd, 




(T • di 
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D 




d,D 











dt 


l^ 
















(J 









-dj + X 


2(7 X 



In components, this Hamiltonian matrix expression becomes. 





mi 




mjdi + djmi 


Dd, 






d 


D 




d,D 










di 


7." 




7/9,+7,", 












aa 




djCTa 





~5^d, 


+ 



crpdi 




Sh/Smj 
Sh/SD 

Sh/Scr 



Sh/Smj 
Sh/SD 
Sh/Sjf 
Sh/Sap 



(2.81) 



(2.82) 



where t"^^ = 2eai3K represents (twice) the vector cross product for liquid crystals. In this expression, 
the operators act to the right on all terms in a product by the chain rule and, as usual, the 
summation convention is enforced on repeated indices. At this point we have switched to using 
both lower and upper Greek indices for the internal degrees of freedom, so that we will agree later 
with the more general theory, in which upper Greek indices refer to a basis set in a Lie algebra, and 
lower Greek indices refer to the corresponding dual basis. Lower Latin indices still denote spatial 
components. 



Remarks about the Hamiltonian matrix. The Hamiltonian matrix in equation ( p. 82 ) was 
discovered some time ago in the context of spin-glasses and Yang-Mills magnetohydro dynamics 
(YM-MHD) by using the Hamiltonian approach in Holm and Kupershmidt [1988]. There, it was 
shown to be a valid Hamiltonian matrix by associating its Poisson bracket as defined in equation 
(|2^ ) with the dual space of a certain Lie algebra of semidirect-product type that has a generalized 
two-cocycle on it. The mathematical discussion of this Lie algebra and its generalized two-cocycle 
is given in Holm and Kupershmidt [1988]. A related Poisson bracket for spin glass fluids is given 
in Volovik and Dotsenko [1980]. A different Poisson bracket for nematic liquid crystals is given in 
Kats and Lebedev [1994], who discuss a constrained Poisson bracket that in general does not satisfy 
the Jacobi identity. The liquid crystal equations in Kats and Lebedev [1994] also differ from the 
Ericksen-Leslie equations by being first order in the time derivatives of the director, rather than 
second order, as in the Ericksen-Leslie theory. See also Isaev et al. [1995] for a discussion of this 
first order theory using the Poisson bracket approach. The present work ignores the first order 
(kinematic) theory in what follows and concentrates on second order (dynamic) theory. 
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Being dual to a Lie algebra, our matrix in equation ( p.82| ) is in fact a Lie-Poisson Hamil- 
tonian matrix. See, e.g., Marsden and Ratiu [1999] and references therein for more discussions 
of such Hamiltonian matrices. For our present purposes, its rediscovery in the PCF context links 
the physical and mathematical interpretations of the variables in the theory of PCFs with earlier 
work in the gauge theory approach to condensed matter, see, e.g., Kleinert [1989]. These gauge 
theory aspects emerge upon rewriting the Lie-Poisson Hamiltonian equations in terms of covariant 
derivatives with respect to the space-time connection one- form given by vdt + 'y^dxm, as done 
in Holm and Kupershmidt [1988]. The gauge theory approach to liquid crystal physics is reviewed 
in, e.g., Trebin [1982], Kleman [1983], and Kleman [1989]. 

The generalized two-cocycle in the Hamiltonian matrix ( |2.82| ) is somewhat exotic for a classical 



fluid. This generalized two-cocycle consists of the partial derivatives in equation ( 2.82 ) appearing 



with the Kronecker deltas in the cr — 7 cross terms. The first hint of these terms comes from 



the exterior derivative dv appearing in the kinematic equation ( |2.62| ) for 7^. Finding such a 
feature in the continuum theory of liquid crystals may provide a bridge for transfering ideas and 
technology between the classical and quantum fluid theories.^ See Volovick [1992] and Volovick 
and Vachaspati [1996] for discussions of similar opportunities for technology transfer in the theory 
of superfluid Helium. The implications of the generalized two-cocycle for the solutions of the liquid 
crystal equations can be seen by considering two special cases: static solutions with one-dimensional 
spatial dependence; and spatially homogeneous, but time-dependent, liquid crystal dynamics. 

Static perfect liquid crystal solutions with z-variation. Static (steady, zero-velocity) solu- 
tions, with one-dimensional spatial variations in, say, the z-direction obey equations ( |2.75| ) - (|2.78 ) 
specialized to 

d 5h d 5h d 5h 

dz 5D ^ dz 573 dz 5a 

d 6h 6h 6h , . 

— — = 273X— +2crx— , 2.83 
dz 073 073 da 

d 6h 5h 
-r ^ = 2 73 X -— . 
dz da oa 



The sum of terms in the first equation of the set ( 2.83| ) vanishes to give zero pressure gradient 



as a consequence of the latter two equations. We compare the latter two equations with the E(3) 
Lie-Poisson Hamiltonian systems, given by 

dn dH ^ dH ^ , 
dV dH ^ 

= — X r . 

dt dn 

When the Hamiltonian H(n, T) in these equations is specialized to 

H=(i:^^)+M9X-r, (2.85) 

^ i=i * ^ 



*In quantum field theory, these partial derivative operators in the Poisson bracket (or commutator relations) are 
called non-ultralocal terms, or Schwinger terms, after Schwinger [1951, 1959]. These terms lead to the so-c alled 



"quantum anomalies." Of course, no quantum effects are considered here. However, the Poisson bracket (2.82) still 
contains classical Schwinger terms. 
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where /j, i = 1,2,3, and Mgx are constants, then the E{3) Lie-Poisson equations (2.84) speciaUze 
to the classical Euler-Poisson equations for a heavy top, discussed in, e.g., Marsden and Ratiu 
[1999]. 

By comparing these two equation sets, one observes that (at least for algebraic h) the static 



perfect liquid crystal flows with z- variation in equations ( 2.83 ) and the E{3) Lie-Poisson tops 
governed by equations ( 2.84| ) are Legendre duals to each other under the map, 



d_ 
dz 



d_ 
di 



/i(73,'T) = n.73 + r-tT-/7(n,r) 



so that 



dh 

Hence, we arrive at the result: 
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dh 
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dH 
dU 
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(2.^ 



(2.87) 



The class of static one dimensional flows of a perfect liquid crystal is Legendre-isomorphic 
to the class of E(3) Lie-Poisson tops. 

These tops conserve 11 • T and |r|^, but in general they are not integrable. 



Spatially homogeneous, time-dependent perfect liquid crystal flows. 

neous solutions of equations (2.75) - ( 2.7^ ) obey the dynamical equations. 



1 da 

2 a 
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6h 



6h 
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6h 



X Ir. 
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Spatially homoge- 

(2.88) 
(2.89) 



For a single value of the spatial index, say m = 3, these are nothing more than the E{3) top 
equations ( 2.84| ) with time re-parameterized by i — > —2t. Hence, in this case, the Hamiltonian 
internal dynamics of a spatially homogeneous liquid crystal is essentially identical to the E{3) Lie- 
Poisson dynamics of a top. In the multi-component case, one sums over m = 1, 2, 3, in the second 
term of equation ( ^.88 ) and, thus, the resulting dynamics is more complex than the simple top. 
Hence, we have: 

The class of spatially homogeneous, time- dependent perfect liquid crystal flows is iso- 
morphic to the generalization ( 2.83^ - (2.8§) of the E{3) Lie-Poisson tops. 



Action principle ^4 — Clebsch representation. Another representation of Hamilton's prin- 
ciple for liquid crystals in the Eulerian fluid description may be found by constraining the Eulerian 
action S in equation (2.55) by using Lagrange multipliers to enforce the kinematic equations (2.61) 



and ( 2.62 ). The constrained action is, thus, 



j dt j d^x i{u,D,v,y) + 



\ dt 



+ 



+ /3, 



' / rn du ^ (97„ duk \ 



dxj 



(2.90) 
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with Lagrange multipliers 4> and (3,^. Stationarity of S under variations in and u implies the 
relations 

^v. ^ + pi + 2/3^X7^ = 0. (2.92) 

These are the Clebsch relations for the momentum of the motion = 6i/6uk and the director 
angular momentum of the micromotion a = bijbv . Stationary variations of 5 in D and 7„ give, 
respectively, the dynamical equations for the canonical momenta, tto = (j) and tt^^ = as 



Finally, variations in the Lagrange multipliers ip B^nd imply the kinematic equations ( 2.61| ) and 



( 



2.62), respectively. These two kinematic equations combine with the four variational equations 
I91I ) through ( p^ ) o recover the Eulerian motion and micromotion equations after a calculation 
using the Clebsch relations, the Kelvin-Noether form of the motion equation, and the dynamical 
equations for the Clebsch potentials. 

At the end of in Section ^, we shall systematize this type of calculation and, thus, clarify its 
meaning as a Poisson map. For now, we simply remark that the evolutionary Clebsch relations are 



Hamilton's canonical equations for the Hamiltonian obtained from the constrained action in (2.90) 



by the usual Legendre transformation. Perhaps not unexpectedly, this Hamiltonian agrees exactly 
with that in equation ( |2.73| ) obtained from the Legendre transformation in u and u alone. 

As we shall discuss more generally in the next section, the Clebsch representations for the mo- 
mentum and director angular momentum provide a Poisson map from the canonical Poisson bracket 
in the Clebsch variables to the Lie-Poisson bracket for the Hamiltonian matrix with generalized 
two-cocycle found in equation (|2.82| ). 

Historically, the Hamiltonian approach has been very fruitful in modeling the hydrodynamics 
of complex fluids and quantum liquids, including superfluids, going back to the seminal work of 
Khalatnikov and Lebedev [1978], Khalatnikov and Lebedev [1980] and Dzyaloshinskii and Volovick 
[1980]. The Clebsch approach has provided a series of physical examples of Lie-Poisson brackets: for 
superfluids in Holm and Kupershmidt [1982]; superconductors in Holm and Kupershmidt [1983a]; 
Yang- Mills plasmas (chromohydrodynamics) in Gibbons, Holm and Kupershmidt [1982], and Gib- 
bons, Holm and Kupershmidt [1983]; magnetohydrodynamics, multifluid plasmas, and elasticity, 
in Holm and Kupershmidt [1983b]; Yang- Mills magnetohydrodynamics in Holm and Kupershmidt 
[1984]; and its relation to superfluid plasmas in Holm and Kupershmidt [1987] and spin-glasses in 
Holm and Kupershmidt [1988]. Many, but not all, of these Lie-Poisson brackets fit into the present 
Euler-Poincare framework for PCFs. 

The Euler-Poincare framework also accomodates many of the various types of Poisson brack- 
ets (such as "rigid body fluids") studied over the years by Grmela, Edwards, Beris, and others, 
as summarized in Beris and Edwards [1994]. For liquid crystals, these authors develop a bracket 
description both for Ericksen-Leslie equations and the Doi-Edwards theory based on the confor- 
mation tensor C, which is related to the director theory by C = n n. The extension of the 
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present results to this case may be accomplished, e.g., by following the Clebsch approach of Holm 
and Kupershmidt [1983b] who treated the corresponding case of Lie-Poisson brackets for nonlin- 
ear elasticity. The treatment in Beris and Edwards [1994] ignores the geometrical content of the 
Lie-Poisson formulation in preference for its tensor properties alone. 



2.4 Summary for perfect liquid crystals 

We now recapitulate the steps in the procedure we have followed in deriving the Euler-Lagrange- 
Poincare-Clebsch equations and the Lie-Poisson Hamiltonian formulations of the dynamics of per- 
fect liquid crystals. 

1. Define the order parameter group and its coset space. 

2. Write Hamilton's principle in the Lagrangian fluid description. 

3. Make 2 stages of reduction: 

1**, to introduce the reduced set of variables in the Lagrangian fluid description; and 
2"^'^, to pass to the Eulerian fluid description. 

4. Legendre transform to obtain the Hamiltonian formulation. 

The alternative Clebsch procedure starts directly with an action for Hamilton's principle that is 
defined in the Eulerian fluid description and constrained by the Eulerian kinematic equations. Its 
Hamiltonian formulation is canonical and passes to a Lie-Poisson formulation via the Poisson map 
that is defined by the Clebsch representations of the momentum and internal angular momentum 
in equations ( p.91| ) and ( p. 92 ), respectively. 



Many physical extensions of these results for perfect liquid crystals are available, e.g., to in- 
clude MHD, compressibility, anisotropic dielectric and diamagnetic effects, linear wave excitation 
properties, etc. However, we wish to spend the most of the rest of this paper setting the formula- 
tions we have established here for perfect liquid crystal dynamics into the geometrical framework of 
Lagrangian reduction by stages developed in Cendra, Marsden and Ratiu [1999]. This geometrical 
setting will take advantage of the unifying interpretation of order parameters as coset spaces of bro- 
ken symmetry groups. (The coset interpretation of order parameters for liquid crystals, superfluids 
and spin glasses is reviewed, e.g., in Mermin [1979].) The present formulations are geometrical 
variants of the Ericksen-Leslie equations for liquid crystal dynamics that illuminate some of their 
mathematical features from the viewpoint of Lagrangian reduction. 



3 Action principles and Lagrangian reduction 

As we have seen, the passage to reduced variables v and 7^ for liquid crystals restricts the variables 
n, n, Vn to the coset space SO{3)/0{2) of rotations that properly affect the director n and imposes 
invariance of the theory under the reflections n — n. The reduced variables transform properly 
under SO{3), because the 0(2) isotropy subgroup of n has been factored out of them. Thus, we 
"mod out" or "reduce" the symmetry-associated degrees of freedom by passing to variables that 
transform properly under rotations in SO{3) and admit the Z2 reflections n — > — n. The removal 
of degrees of freedom associated with symmetries is the essential idea behind Marsden- Weinstein 
group reduction in Marsden and Weinstein [1974]. Marsden- Weinstein reduction first appeared 



Dynamics of perfect complex fluids 



August 1, 2000 



24 



in the Hamiltonian setting. However, this sort of reduction by symmetry groups has been recently 
extended to the Lagrangian setting, see Cendra, Marsden and Ratiu [1999] and Marsden, Ratiu and 
Scheurle [1999]. The remainder of the paper apphes the mathematical framework of Lagrangian 
reduction by stages due to Cendra, Marsden and Ratiu [1999] to express some of the properties 
of PCF dynamics in the Eulerian description for an arbitrary order parameter group. 

A synthesis of the nonlinear dynamics for the motion and micromotion of various perfect complex 
fluid models is possible, due to their common mathematical basis. The mathematical basis common 
to all ideal fluid motion - both classical and complex fluids - is Hamilton's principle, see, e.g., 
Serrin [1959], 

6S = 6 j Ldt = 0. (3.1) 

In the Lagrangian (or material) representation for fluids, the motion is described by the Euler- 
Lagrange equations for this action principle. 

In the Eulerian (or spatial) representation for fluids, the Euler-Lagrange equations for the 
dynamics are replaced by the Euler-Poincare equation. The distinction between Euler-Lagrange 
equations and Euler-Poincare equations is exemplified by the distinction between rigid body motion 
expressed in terms of the Euler angles and their time derivatives on the tangent space TSO{3) of 
the Lie group of proper rotations 50(3), and that same motion expressed in body angular velocity 
variables in its Lie algebra so{3). Poincare [1901] was the first to write the latter equations on an 
arbitrary Lie algebra; hence, the name Euler-Poincare equations. 

Euler-Poincare equations may be understood and derived via the theory of Lagrangian reduction 
as in Cendra et al. [1999], and Cendra, Marsden and Ratiu [1999]. Euler-Poincare equations arise 
when Euler-Lagrange equations and their corresponding Hamilton principles are mapped from a 
velocity phase space TQ to the quotient space TQ/G (a vector bundle) by a Lie-group action 
of a symmetry group G on the configuration space Q. If L is a G- invariant Lagrangian on TQ, 
this process maps it to a reduced Lagrangian and a corresponding reduced variational principle 
for the Euler-Poincare dynamics on TQ/G in which the variations are constrained. See Weinstein 
[1996] and Cendra, Marsden and Ratiu [1999] , for expositions of the mathematical framework that 
underlies Lagrangian reduction by stages and Holm, Marsden and Ratiu [1998] for a discussion 
of Euler-Poincare equations and their many applications in classical ideal fluid dynamics from the 
viewpoint of the present paper. See Marsden, Ratiu and Scheurle [1999] for additional insight and 
recent results in Lagrangian reduction. 

The order parameters of PCFs are material variables. The Lagrangian in Hamilton's principle 
(|3T|) for PCFs is the map, 

L : TG X y* X TO I — >R. (3.2) 

That is, the velocity phase space for the PCF Lagrangian L in material variables is the Cartesian 
product of three spaces: 

TG, the tangent space of the Lie group G of fluid motions (the diffeomorphisms that take the fluid 
parcels from their reference configuration to their current positions in the domain of flow), 

V*, the vector space of advected quantities carried with the fluid motion, and 
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TO, the tangent space of the Lie group O of fluid micromotions {O is the order parameter Lie 
group) . 

The advected quantities in V* include the volume element or mass density and whatever else is 
carried along with the fluid parcels, such as the magnetic field intensity in the case of magneto- 
hydrodynamics. The new feature of PCFs relative to the simple fluids with advected parameters 
treated in Holm, Marsden and Ratiu [1998] is the dependence of their Lagrangian on TO. The 
order parameter coset space at each material point is acted upon by the order parameter Lie group. 
(We choose the convention of group action from the right.) Since the order parameter is a ma- 
terial property, the diffeomorphism group G also acts on the order parameter group, as 
O X G — > O, denoting action from the right. 

In this Section, we shall use Hamilton's principle (3.1) with Lagrangian (3.2) to obtain the 



dynamical equations for the motion and micromotion of PCFs whose order parameters are defined 
as coset spaces of Lie groups. In doing so, we shall begin by assuming this Lagrangian is invariant 
under the right action of the order parameter Lie group O on its tangent space TO. (This right 
action on the space of internal variables leaves the other components of the configuration space 
TG and V* fixed.) We shall assume this Lagrangian is also invariant under the right action of 
the diffeomorphisms G, which relabel the fluid parcels. (This action of G does indeed affect the 
material variables defined on TO and V*.) Under these symmetry assumptions we shall perform 
the following two group reductions 

(TG xV* X {TO/0))/G ~ g X (y* X o)g-^{t) , 

with respect to the right actions of first O and then G, by applying group reduction to the velocity 
phase space of this Lagrangian in two stages, 

1st stage, (TG xV* x T0)/0 - TG x V* x o , (3.3) 
2nd stage, (TG x y* x o)/G ~ g x {V* x o)g~^{t) . (3.4) 

Here we denote isomorphisms as, e.g., o ~ TO/0 and g ~ TG/G, where the Lie algebras o and 
g correspond, respectively, to the Lie groups O and G. The first stage is Lagrangian reduction by 
the right action of O, the order parameter group |^ The second stage is Lagrangian reduction of 
the first result by the right action of the diffeomorphisms G in the first factor and by composition 
of functions in the second factor. Because of the assumed invariances of our Lagrangian, these 
two stages of reduction of the velocity phase spaces will each yield a reduced Lagrangian and a 
corresponding reduced variational principle for the dynamics. The group actions at each stage are 
assumed to be free and proper, so the reduced spaces will be local principle fiber bundles.^ The 
mathematical formulation of the process of Lagrangian reduction by stages and the introduction 
of various connections on the Lagrange- Poincare bundles that arise in Lagrangian reduction 
are discussed in Cendra, Marsden and Ratiu [1999]. These Lagrange-Poincare bundles are special 
cases of Lie algebroids. See Weinstein [1996] for a fundamental description of the relation between 
Lagrangian mechanics and Lie algebroids. 

■''Variants exist. In particular, for liquid crystals, only a part of the Lie algebra o is required; namely, so(3)/0(2), 
the part of the Lie algebra so(3) that is invariant under the 0(2) isotropy group of the director, n. 

®A natural fiat connection appears on this bundle, but this bundle picture should be made intrinsic and global, 
while including defect dynamics. A strategy for obtaining equations for the defect dynamics will be discussed in the 
next Section. However, the global bundle picture is for future work. 
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3.1 Lagrangian reduction by stages 

We are dealing with a Lagrangian defined by the map 

L{g,g,ao,X,dx):TGxV* xTO^R, 

where G is the difFeomorphism group that acts on both the vector space V* of advected material 
quantities and the order parameter group O. We assume that L has the following invariance 
properties, 

L{9: 9, ao, X, dx) = L{g, 9, clq, x^, dxip) = L{gh, gh, oq/i, x^/i, dxi^h) , (3.5) 

for all "0 € O and h £ G. In particular, we shall choose ^p = x^^(^) in the first stage and h = g~^{t) 
in the second stage of the reduction, so that 

L{9,g,ao,X,dx) = L{9,g,ao,XX~'^,dxX~^) , (3-6) 
after the first stage of reduction, and 

L{9,9,ao,X,dx) = L{e,gg~'^,a,xx~^9~^,dxx^^9'^) = l{^,a,^,l) , (3.7) 

after the second stage, with ^ = gg~^, a = aog^^, v = {xx ^^)9^^ and 7 • = {dxx~'^)9~^ ■ After 
the first stage of reduction, the reduced action principle yields the Lagrange-Poincare equations, 
and after the second stage we shall obtain the Euler-Poincare equations for a perfect complex fluid 
with an arbitrary order parameter group. 

3.1.1 Lagrange-Poincare equations 

The first stage 

TG xV* xTO < — >TGxV* xo, 



of the two-stage symmetry reduction in ( |3.3D - (3^) affects only the internal variables and passes 



from coordinates on the order parameter Lie group, O, to coordinates on its Lie algebra, o, obtained 
from the tangent vectors of the order parameter Lie group at the identity by the isomorphism 
o ~ TO/0. The results at this first stage consist of Euler-Lagrange equations for the fluid motion, 
coupled through additional components of the stress tensor to equations of a type called Lagrange- 
Poincare equations in Marsden and Scheurle [1995], Cendra, Marsden and Ratiu [1999]. In our 
case, these Lagrange-Poincare equations describe the micromotion in the Lagrangian (or material) 
fluid description. 

The first stage of reduction results in the Lagrange-Poincare action principle, 

djdtjd^X L{x,J,i^,j) =0, (3.8) 

written in the material representation and denoted as follows, 
L{x,J, u, 7) is the reduced Lagrangian on TG x V* x o, 
J{X,t) = det{dx/dX) G V* is the volume element. 
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^{Xjt) = XX ^i^^t) G o is the material angular frequency and 
7 • = dxx^^i^^ t) £ o with components denoted as 7m given by 

7mdXm{X,t) = dxx'^{X,t), 

= lm^dXA=lT\X,t)dXA, (3.9) 

oXa 

where 7 • is the Cosserat strain one-form introduced in Cosserat [1909] and superposed 
"dot" ( )' denotes time derivative at fixed material position X. 

These material quantities satisfy auxiliary kinematic equations, obtained by differentiating their 
definitions, 

{J-^d^xy = {d^X)' = {), (3.10) 
(7 • dx)' = {dxx^^y = dv + a.dv{'y ■ dx) , (3.11) 

The material angular frequency v = xx^^ ^-nd the material Cosserat strain one-form j-dx = dxx~^ 
take their values in the right-invariant Lie algebra o of the order parameter Lie group O. The ad- 
operation appearing in equation ( 3.11 ) denotes multiplication, or commutator, in the Lie algebra 
o. 

The dynamical Lagrange- Poincare equations determine the complex fluid's motion with 
fluid trajectory (^t(X) = x{X,t) with 4>t £ G and micromotion xi^^t) £ O \n the material fluid 
description. These equations take the following forms, 

j-ifdL\^ d dL d / j-i dL ^\ _ ^ ^2) 
V dip J dxp dJ dxm \ d^m ' ^ / ' 

f^V _ ad*— + J^fj-i^V ad* ^ = (313) 

The ad* -operation appearing in ( |3.13| ) is defined in terms of the ad-operation and the symmetric 
pairing (J) : o* x o ^ R between elements of the right Lie algebra o and its dual o* as, e.g.. 

In a Lie algebra basis satisfying [eQ,,e^] = t^^e^ and its dual basis e'^ satisfying (e'',et) = 6^^, we 
may write this formula as 



ad:^ , S ) = ( , ad.S ) = ^ z." (3.15) 



dL \ /dL \ dL 

-d^'^/ = \-d^''''^^^l = d^^'-P' 

Thus, for the sign conventions we choose in (3.14), the ad*-operation is defined as the negative 
transpose of the ad-operation. 

The dynamical equations ( 3.12| ) and ( |3.13| ) follow from the Lagrange-Poincare action principle 
for PCF dynamics in the material fluid description. These Lagrange-Poincare equations may 
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be calculated directly, as 

= sfdtfd^X L{x,J,i^,j) 



- M 


d^X 


- H 


d^X 1 




:-( 




(>d^S 



dL ^. dL^^ /dL ^ . 
-^OXp + -^dJ + ( — ,du } + 
OXp oJ \oiJ 



dL 

dXr. 



J 9 dL ^ 



dL 

d 



j +ad 
dXA 



j- 



dxp dJ 



d 



dXr 



V dx 
dL 



j^i dL 
dL 



d^m '' 



dL \ ^ dL ^ 

,7 )dxp + J—dXr 



d% 



dJ 



Here we define Ti = 5xx ^, in terms of which we calculate, 

5v = Ti — \v, 

dy 



S — adjy S , 



Since 7™ = 7^"* (5X^/5, 



dXA 
this means that 
dXB d 



dT. 
dXA 



cld^ mat X] , 

'A 



-7p 



dXrn dXs " ^ 

dXB d ^ , 
dxm dXs 



Sx^ + ^6jr' 



dx 

dXAdT^_ , ^ 
dx^dXA 



(3.16) 



(3.17) 
(3.18) 



(3.19) 



Thus, the variations in 7^ couple the two Lagrange-Poincare equations. We also drop endpoint 
terms that arise from integrating by parts in time, upon taking 5xp and S to vanish at these 
endpoints. The natural boundary conditions 



dC^ 
dJ 



and fin 



dC 

dim 



0. 



(3.20) 



ensure that the fluid pressure and the normal stress are continuous across a fluid interface. 



3.1.2 Euler-Poincare equations 

The passage next from the Lagrangian fluid description of continuum mechanics to the Eulerian 
fluid description will yield the Euler-Poincare equations. We obtain these equations by applying to 
Hamilton's principle the second stage, 

TG X y* X o I — . g X (y* X o)g-^{t) , 

of the two-stage Lagrangian reduction in ( |3.3| ) - ( |3.4D . This second stage of reduction results in the 
Euler-Poincare action principle, 



6 / /(^, a, i^, 7) = , 



(3.21) 
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with constrained variations 

= 5a = -ar], 

dYl 

Si' = — VS-ad^S-i^r?, (3.22) 
S (7 ■ dx) = dT, + ads (7 ■ dx) — (7 • dx) ij , 

where r](t) = 6g(t)g{t)~^ G g and = (^x(^)x(^)^^ € o both vanish at the endpoints. 

The Euler-Poincare action principle produces the fohowing equations defined on g x {V* x 
o)g~^{t) for the motion and micromotion, in which d/dt denotes Eulerian time derivative at fixed 
spatial position x, 



d SI SI SI SI SI 

dtT^ = -"^«^ + ^^"+;^""+^^^-' 

d SI ,. /^Sl Sl\ SI SI 

These are the Euler-Poincare equations for a perfect complex fluid. In these equations, I is the 
reduced Lagrangian on g x {V* x o)g~^{t). Also, 

ad5»7=[C,??], with C,?7Gg, (3.25) 

is the commutator in the Lie algebra of vector fields, g. In addition, we define the two operations 
ad| and o as 



ad||,r?\ ^ -/|,ad,,?), (3.26) 



and 



The concatenation ar] denotes the right Lie algebra action of G g on a G (by Lie derivative). 
The pairing (• , •) now includes spatial integration and, thus, allows for integration by parts. Similar 
definitions hold for {SI / S^m o 7m) and {Sl/Suov). 
In components, these quantities are given by, 



SI \ SI \ d ( SI 3\ SI dy, 



^^Im I \Slm I \(^Xm\Sj^ V Sjl^ 

/ SI \ I SI \ I SI dv^ \ 



' SI ^ \ I SI ^ \ I ^ d SI ^ 
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Remark. At this point, one might have also introduced (3 + 1) covariant derivatives acting 
on o— valued functions of space and time, as 

d d 
Dm = -Q^ ad^m > and Dt = - ad^ , (3.29) 

with associated curvature (or Yang-Mills magnetic field) given by 



ads^ 



(3.30) 



whose components are expressed, cf. equation (|2.2^ ), 

5i" = <,-7i:^ + t;?«7f7/, (3.31) 

as in Holm and Kupershmidt [1988]. However, the operations ad, ad*, o and Lie derivative are 
sufficient for our present purposes. 

Eulerian kinematic equations. By definition, an Eulerian advected quantity a € V*g{t)~^ 
satisfies 

da 

This advection relation may be written equivalently as 

- + i:,a = o, 

where £^ is the Lie derivative with respect to ^ = g{t)g~^{t), the Eulerian fluid velocity, often 
denoted also as u{x,t). 

The Eulerian versions of the Lagrangian kinematic equations ( 3.1C| ) and ( 3.11 ) are given in 
terms of the Lie derivative by 

— + £^^(^Dd^x^ = 0, (3.32) 
+ (7 • dx) = du + ad,y(7 • dx) . (3.33) 

In these equations, the quantity D{x,t) = J^^{X,t)g{t)^^ is the Eulerian mass density and the 
quantities ij{x,t) = i'{X,t)g{t)^^ and 

7(x,t) • dx = {j{X,t) ■ ^{X,t) ■ dx)g{t)-^ 

are the Eulerian counterparts of the right-invariant material quantities ^{X, t) and 7(X, t) in equa- 
tions (|1|). 
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Remark. We note that equation ( p.33|) implies the — circulation theorem for PCFs, cf. equa- 
tion ( ^l63|) , 



It 



7 • dx 



c{0 



ad,y(7 • dx) . 



(3.34) 



Thus, the circulation of 7 around a loop c(^) moving with the fluid is conserved when adi, 7 is a 
gradient. Otherwise, the curl of this quantity generates circulation of 7 around fluid loops. 

The Euler-Poincare equations ( ^.23 ) and ( 3.24 ) may be obtained directly from the Euler- 
Poincare action principle (3.21), as follows. 



Euler-Poincare action variations. We compute the variation of the action ( ^.21 ) in Eulerian 
variables at fixed time t and spatial position x as. 




5a 



51 



— — - adf 7- + ao — + 1/0 — +7mO 
5^ 



5l_ 
5^ 

5J_ 

5a 



5u) + 

51 
5^ 



51 

5lm 



51_ _51_ 



51 



dtld^xl- 



51 



51 



+ 



d 



dXr 



Vj 6 



5uf^ 
51 51 



+ ad*-^ + ad* ^ 
51^ ^™ 5% 




(3.35) 



56 



5^i 



+ 



51 51 



5u(^ 



5jm 



where we have used the variational expressions in ( ^.22 ) and integrated by parts. Here 77 J a in the 
first boundary term denotes substitution of the vector field r] into the tensor differential form a and 
(• , •) denotes the natural pairing between a tensor field and its dual. 

The dynamical Euler-Poincare equations ( ^.23 ) and ( ^.24 ) are thus obtained from the Euler- 
Poincare action principle ( |3.21 ), by requiring the coefficients of the arbitrary variations 77 and S to 
vanish in the variational formula ( 3.35| ). The remaining terms in ( [3.35 ) yield Noether's theorem 
for this system, which assigns a conservation law to each symmetry of the Euler-Poincare variational 
principle. 



Momentum conservation. In momentum conservation form, the PCF motion equation in the 
Eulerian fluid description ( |3.23| ) becomes, for algebraic dependence of the Lagrangian density / on 

{(,a,U,-frn), 

d dl d f ^ dl ,^ dl p\ J d , dl 



In this equation, expressed in Cartesian coordinates, there is an implied sum over the various 
types of advected tensor quantities, a. This momentum conservation law also arises from Noether's 
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theorem, as a consequence of the invariance of the variational principle ( p.21| ) under spatial trans- 
lations. In fact, the simplest derivation of this equation is obtained by evaluating the variational 
formula ( |3.35| ) on the equations of motion and using the translational symmetry of the Lagrangian 
in Noether's theorem with rjj = d/dxj. For an algebraic Lagrangian density -D, z^, 7m), this 



momentum conservation law becomes, cf. equation (|2.6 



d_dl_ 



d 



dXr. 



dl 



dl 



p '3 



(3.37) 



Kelvin-Noether circulation theorem for PCFs. Rearranging the motion equation ( |3.23| ) and 
using the continuity equation for D in (p.32|) gives 



di 



1 61 



dxA 



c{0 



1 

15 



61 



61 61 
— oa + — ov + — 
da dl' 07. 



(3.38) 



where the circulation loop c{(,) moves with the fluid velocity ^ and we have used the following 
relation, valid for one-form densities, 



ad2 



61 



£ 



61 



6^ -~^'6r ys^j 

This relation may be checked explicitly in Cartesian coordinates, as follows. 



with 



61 



61 



dxj (g) d X I . 



(3.39) 
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dx^V 6^^ 
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61 9gJ 
W'dx'' 



rf d^x 



6C^ V d^' dx'' 



(|,ad^ry 



ad^ - , 7? 



(3.40) 



See Holm, Marsden and Ratiu [1998] for more explanation and discussion of the Kelvin-Noether 
circulation theorem for Euler-Poincare systems. 

In components, the Kelvin-Noether circulation theorem ( |3.38 ) for PCFs with Lagrangian /(^, D, u, 
may be written using equation ( |3.28| ) as. 



dt 



HO 



1 61 , 

dxi 

D6i, ^ 



Hi) 



D 



D 



d 61 
dxj 6D 



61 dv^ 
dul^ dx- 



+ 



d 



dXr. 



61 



67!, 



51 97^ 



6% 



I dxi 



dxj. 



(3.41) 



Thus, gradients of angular frequency and order parameter strain may cause fluid circulation. 



3.2 Hamiltonian dynamics of PCFs 

The Legendre Transformation. One passes from Euler-Poincare equations on a Lie algebra g 
to Lie-Poisson equations on the dual g* by means of the Legendre transformation, see, e.g., Holm, 
Marsden and Ratiu [1998]. In our case, we start with the reduced Lagrangian Z on g x {V* x o)g{t)~^ 
and perform a Legendre transformation in the variables ^ and v only, by writing 

61 61 

/x = — , cr = — , h{fi,a,a,-f) = {fi,C) + -l{^,a,u,-f). (3.42) 



5^ 



6v 
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One then computes the variational derivatives of h as 



6h 
6fj, 



dh 
6^ 



Sh 
5a 



6l_ 6h 

6a ' (57 



6l_ 

6j 



(3.43) 



Consequently, the Euler-Poincare equations ( 3.23 ) - ( |3.33| ) for PCF dynamics in the Eulerian de- 
scription imply the following equations, for the Legendre-transformed variables, (/x,a, cr, 7), cf. 
equations (|77|) - (|]78|) 



97 
'dt 



'dt 
da 

at 

• dx 

da 
'dt 



5h Sh 
ad^;,/^^ /i - ^ o a - — X. 7m 



6h 
6^ 



o a . 



£ 



£5h/Sf,il ■ dx) + d(^—^ + adsh/Sail ■ dx) , 



(3.44) 



As for the case of liquid crystals discussed earier, these equations are Hamiltonian and may be 
expressed in terms of a Lie-Poisson bracket. 



Lie-Poisson bracket for PCFs. Assembling the PCF equations ( |3.44| ) into Hamiltonian form 
gives, symbolically. 
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dt 
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— (grad — ad-y)n 
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-(div-ad;)n 
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5h/ 5^ 




Sh/da 




Sh/Sj 




Sh/Sa 



(3.45) 



with boxes □ indicating where the matrix operations occur. In the 7 — cr entry of the Hamiltonian 
matrix (3.45), one recognizes the covariant spatial derivative defined in equation ( 3.2S| ), and 
finds its adjoint operator in the cr — 7 entry. More explicitly, in terms of indices and differential 
operators, and for a = D, the mass density, this Hamiltonian matrix form becomes 





Mi 








dnt - 


7^- 




d 


D 

















dt 


It 




7fa. + 7,", 























-5id, + 




+ K „ 



6h/5fij 
6h/SD 
Sh/Sjf 
Sh/Sap 



(3.46) 



Here, the summation convention is enforced on repeated indices. Upper Greek indices refer to 
the Lie algebraic basis set, lower Greek indices refer to the dual basis and Latin indices refer to 
the spatial reference frame. The partial derivative dj = d/dxj, say, acts to the right on all terms 
in a product by the chain rule. For the case that are structure constants e^^re for the Lie 
algebra so(3), the Lie-Poisson Hamiltonian matrix ( |2.82| ) for liquid crystals is recovered, modulo 
an inessential factor of 2. 



Remark. As mentioned earlier in our discussion of Hamiltonian dynamics of liquid crystals, the 
Hamiltonian matrix in equation ( 3.46 ) was discovered some time ago in the context of investigating 
the relation between spin-glasses and Yang-Mills magnetohydrodynamics (YM-MHD) by using 
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the Hamiltonian approach in Holm and Kupershmidt [1988]. There, it was shown to be a vahd 
Hamiltonian matrix by associating its Poisson bracket with the dual space of a certain Lie algebra 
of semidirect-product type that has a generalized two-cocycle on it. This generalized two-cocycle 
contributes the grad and div terms appearing in the more symbolic expression of this Hamiltonian 
matrix in equation ( |3.45| ). 

The mathematical discussion of this Lie algebra and its generalized two-cocycle, as well as the 
corresponding Lie- Poisson Hamiltonian equations for spin-glass fluids and YM-MHD, are given in 
Holm and Kupershmidt [1988]. The present work provides a rationale for the derivation of such 
Lie-Poisson brackets from the Lagrangian side. 



Spatially one-dimensional static solutions with z-variation. Static (steady, zero-velocity) 
solutions for PCFs, with constant pressure and one-dimensional spatial variations in, say, the z- 
direction obey equations ( p. 46 ) for i = 3 = j, rewritten as 

^ d 5h ^ d 5h d dh ^ 

dz 6D dz 6jf dz 5aa ' 

d 5h ^c, ^ 5h ^ 5h 

d 6h a K 

As for the case of liquid crystals, the sum of terms in the first equation of the set ( 3.47] ) vanishes 
to give zero pressure gradient, as a consequence of the latter two equations. 
Under the Legendre transformation 



/i(73,ct) = n„73" + rV„ - HiU,T) , (3.48) 



these equations become. 



^n. = c||:n, - (3.49) 

Tz^ ~ ^f'^m:^ ■ 

These Legendre-transformed equations are Poincare's [1901] generalization of Euler's equations for 
a heavy top, expressing them on an arbitrary Lie algebra with structure constants t^^. Thus, 

The steady, spatially one- dimensional solutions for all PCFs have the underlying Lie 
algebra structure discovered in Poincare [1901]. 



Spatially homogeneous, time-dependent PCF flows. Spatially homogeneous solutions of 
equations ( |3.46| ) obey the dynamical equations, 

dca PI 6h ^ Sh , , 
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For a single spatial index, say 7^, these are again Poincare's [1901] equations generalizing Euler's 
equations for a heavy top to an arbitrary Lie algebra. Of course, the corresponding Hamiltonian 
matrix for this system is the lower right corner of the matrix in equation (|3.46| ) . 
When f^^ 

— Co/3k for the Lie algebra so(3), Poincare's equations (|3.50[) correspond to the 
Leggett equations for ^He — A with spin density ctq and spin anisotropy vector jf, see Leggett 
[1975]. For special solutions of these and other related equations in the context of ^He — A, see 
Golo and Monastyrskii [1977], Golo and Monastyrskii [1978], Golo et al. [1979]. 

Evolution of the disclination density. Holm and Kupershmidt [1988] use the chain rule and 



the defining relation for the disclination density, cf. equation (2.28) 



(3.51) 



to transform the Hamiltonian matrix ( 3.46| ) to a new Hamiltonian matrix, whose Lie-Poisson Hamil- 
tonian dynamics may be written as 



d_ 

dt 



D 
B 
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(3.52) 
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Hkdi + dkiJLi 
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-^;m,4 + d-niBd 


dkD 








+ Btkdj-B^^d, 








dkCTa 








diBt 





4- a r> K 



Sh/S^k 
Sh/dD 

Sh/SBH 



The corresponding PCF dynamics for the disclination density emerges as 



dt 



6h 



6h 



(3.53) 



As expected, this equation preserves the trivial solution = 0, which is the case when the 
7m— strain field is continuous and the complex fluid has no defects. We shall mention a strategy 
for dealing with defects in Section ^ 



3.3 Clebsch approach for PCF dynamics 

Following Serrin [1959], we call the auxiliary constraints imposed by the Eulerian kinematic equa- 



tions the Lin constraints. As we shall see, the diamond operation o defined in equation (3.27) 
arises naturally in imposing the Lin constraints. Taking variations of the constrained Eulerian 
action, 

S = jdt !^l{(,a,u,j) + (^v,^ + £i:a) + (^P,^ + £a-du-ad,^)'^ (3.54) 
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yields the following PCF Clebsch relations, 



SI 

: — -voa-Po^ = 0, 



6l_ 

61 dv 



5v: — + d/? - ad; /3 = , 



5a : i^fV = 0, (3.55) 

5a dt ^ ^ ' 

da 

5v : — — h a = , 
at 

5l ■■ ^ - ^ - £5 /3 + ad* /? = , 

SP : — — \- £c'y — di' — ad^ 7 = 0. 
ot 

We shall show that these Clebsch relations recover the Euler-Poincare equations ( p.23| ) - ( 3.241 ). (In 
what follows, we shall ignore boundary and endpoint terms that arise from integrating by parts.) 
The diamond operation o is defined by 

{v o a , rj) = — {v , a) = — {v , a 7]) . (3.56) 

This operation is antisymmetric, 

{v o a , rj) = — {ao V , r]) , (3.57) 

as obtained from, 

{v , £r^a) + {£rjV , a) = , or, {v , arj) + {vrj , a) = , (3.58) 

and the symmetry of the pairing (• , •). The diamond operation also satisfies the chain rule under 
the Lie derivative, 

{£^ {v o a) , rj) = {{£^ v) o a , rj) + {v o {£^ a) , rj) . (3.59) 

This property can be verified, as follows, 

{£^ V o a , rj) + {v o £^a , rj) = {v S^r] , a) — {v rj ^ , a) 

= {a , V (ad^ r])) = — {aov , (ad^ r])) 

= {adl{aov),rj) = {£^{aov),rj), (3.60) 



where we have used {v ^ , a rj) + {v £,rj , a) = 0, implied by ( p.58|) , in the first step. Finally, we have 
the useful identity, 

{(3odv,rj) = -{dl3ov,rj), (3.61) 
as obtained from {dv)rj = d{vrj) and 

{I3,d{urj)) + {df3,vrj)={). (3.62) 
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These three properties of the o operation and the PCF Clebsch relations (|3.55|) together imply 



o a + /3 o 7^ 



61 



51 61 
— oa + — ou + — 
da dv 07; 



(3.63) 



+ 



(ad* /3) o 7 + /? o ad,, 7 + (ad* /5) o 



The term in square brackets is seen to vanish, upon pairing it with a vector field, integrating by 
parts and again using the properties of the o operation. This manipulation recovers the PCF 
motion equation ( 3.23 ) as 



d 



\ 61 61 



61 



dt 



since, as we have seen. 



6a 



6i 



6u 



* 51 



61 



(3.64) 



(3.65) 



for one- form densities such as 61/6^. 

The PCF micromotion equation (3.24) is also recovered from the Clebsch relations ( |3.55| ). 
This is accomplished by taking the time derivative of the (^z^— formula, substituting the 6(3— and 
57— formulas, and using linearity of ad* to find 



61 , 61 

h d — 

6v 07 



(3.66) 



(i(ad* /?) + ad^^ (3 - ad* (ad* (5) - ad 

ad* dj3 — ad* (ad* /3) 
61 



= 6. ■ 

Hence, the Clebsch relations ( |3.55 ) also recover the micromotion equation ( |3.24| ). 

Remarks. From the Hamiltonian viewpoint, the pairs (f , a) and (/3, 7) are canonically conjugate 
variables and the Clebsch map (u,a, /3,7) — > {^J-jCr), with 

61 



voa + (3o-f, (3.67) 

— = a = -dp + adZp, (3.68) 
dv ' 

is a Poisson map from the canonical Poisson bracket to the Lie-Poisson Hamiltonian structure 
given in equation (3.46), in which a = D. Of course, there is no obstruction against allowing a to 
be any advected quantity, as discussed in Holm, Marsden and Ratiu [1998]. 

The generalized two-cocycle associated with the Hamiltonian matrix in ( p. 46 ) arises from the 
term dp in the cr— part of this Poisson map. 

Various other applications of the Lin constraint and Clebsch representation approach in for- 
mulating and analyzing ideal fluid and plasma dynamics as Hamiltonian systems appear in Holm 
and Kupershmidt [1983b], Marsden and Weinstein [1983], Zakharov et al. [1985], Zakharov and 
Kusnetsov [1997]. 
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3.4 Conclusions for PCFs 

Perfect complex fluids (PCFs) have internal variables whose micromotion is coupled to the fluid's 
motion. Examples of PCFs include spin-glass fluids, superfluids and liquid crystals. PCF internal 
variables are materially advected order parameters that may be represented equivalently as either 
geometrical objects, or as coset spaces of Lie groups. The new feature of PCFs relative to simple 
fluids with advected parameters treated in Holm, Marsden and Ratiu [1998] is the dependence of 
their Lagrangian 

L:TGxV* xTOi — >IL, 

on TO, the tangent space of their order parameter group. Moreover, the diffeomorphisms G act on 
TO. We treat Lagrangians that are invariant under the right actions of both the order parameter 
group O and the diffeomorphisms G. In this case, reaching the Euler-Poincare fluid description 
requires two stages of Lagrangian reduction, 

{TG xV* X (TO/O)) /G ~ g X (F* X o)g-\t) , 

rather than the single stage of Lagrangian reduction (with respect to the "relabeling transforma- 
tions" of G) employed for simple fluids in Holm, Marsden and Ratiu [1998]. 

After studying the example of nematics in Section ^, we derived the Euler-Poincare dynamics 
of PCFs in two stages of Lagrangian reduction in Section |3[ The first stage produced the Lagrange- 
Poincare equations derived from an action principle defined on the right invariant Lie algebra of 
the order parameter group in the Lagrangian (or material) fiuid description. The second stage of 
Lagrangian reduction passed from the material fluid description to the Eulerian (or spatial) fluid 
description and produced the Euler-Poincare equations for PCFs. We also derived these Euler- 
Poincare equations using the Clebsch approach. 

In addition, we used a Legendre transformation to obtain the Lie-Poisson Hamiltonian formula- 
tion of PCF dynamics in the Eulerian fluid description. The Lie-Poisson Hamiltonian formulation 
of these equations agreed with that found earlier in Holm and Kupershmidt [1987] , Holm and Ku- 
pershmidt [1988], who treated spin-glass fluids, Yang- Mills magnetohydrodynamics and superfluid 
^He and ^He. Thus, we found that Lagrangian reduction by stages provides a rationale for deriving 
this Lie-Poisson Hamiltonian formulation from the Lagrangian side. This approach also fits well 
with some gauge theoretical descriptions of condensed matter physics. 

4 A strategy for introducing defect dynamics 

Many other potential applications of the Euler-Poincare framework abound in the physics of con- 
densed matter. For example, besides the perfect liquid crystal dynamics treated here explicitly, 
the superfluid hydrodynamics of the various phases of ^He may be treated similarly. In particu- 
lar, the geometrical framework of Lagrangian reduction by stages is well-adapted to the standard 
identiflcation of the phases of ^He with the independent cosets of the order parameter group 
50(3) X S0{3) X U{1), as discussed, e.g., in Mineev [1980] and Volovick [1992]. Magnetic mate- 
rials may also be treated this way. The seminal papers on the geometrical properties of magnetic 
materials are Dzyaloshinskii [1977], Volovik and Dotsenko [1980] and Dzyaloshinskii and Volovick 
[1980]. Other recent studies of the dynamics of magnetic materials and superfluid ^He in directions 
relevant to the present paper appear, e.g., in Holm and Kupershmidt [1988], Balatskii [1990], Isayev 
and Peletminsky [1997], Isayev, Kovalevsky and Peletminsky [1997]. 
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Most of these physical apphcations involve defects (imperfections, or "glitches" that appear as 
discontinuities in the order parameter) and one must describe their dynamics, as well. This is an 
area of intense investigation in many contexts in condensed matter physics. Our approach is based 
on an analogy with the theory of the Hall effect in a neutral multicomponent fluid plasma when 
inertia is negligible in one of the components (the electrons). 

The Hamiltonian context for Hall effects in a neutral ion-electron plasma was considered by 
Holm [1987] for a normal-fluid plasma and by Holm and Kupershmidt [1987] for a multicomponent 
electromagnetically charged superfluid plasma. In this context, the introduction of an independent 
gauge fleld associated with the momentum of a distribution of superfluid vortex lines generalizes to 
apply for defects or vortices in any continuous medium possessing an order parameter description 
that arises from spontaneous symmetry breaking. In the remainder of this paper, we shall use the 
idea of reactive forces arising via the Hall effect to discuss the case of quantum vortices in superfluid 
Helium-II. In particular, we shall apply the Hall effect analogy to describe the hydrodynamics of the 
quantum vortex tangle in superfluid turbulence as an additional "third fluid" in Landau's two-fluid 
fluid model. The third fluid associated with the vortex tangle carries momentum and moves with 
its own independent velocity in superfluid flows. 

4.1 Vortices in superfluid ^He 

In superfluid ^He the order parameter group is U{1) and the defects are called vortices. These are 
quantum vortices, since their circulation comes in integer multiples of k = h/m ~ 10~^cm^/sec. 
Conservation of the number of quantum vortices moving through superfluid ^He (and across the 
streamlines of the normal fluid component) is expressed by 

4- / uj-ndS = 0, (4.1) 

dt Js ' ^ ^ 

where the superfluid vorticity u is the areal density of vortices and n is the unit vector normal to 
the surface S whose boundary dS moves with the vortex line velocity V£. When u = curlv^ this is 
equivalent to a vortex Kelvin theorem 

^(f v,-dx = 0, (4.2) 

JdS{ve) 

which in turn implies the fundamental relation 

-^-v^xa; = V/x. (4.3) 

The superfluid velocity naturally splits into = u— A, where u = V0 and (minus) the curl of A 
yields the superfluid vorticity u. The phase (p is then a regular function without singularities. This 
splitting will reveal that the Hamiltonian dynamics of superfluid ^He with vortices may be expressed 
as an invariant subsystem of a larger Hamiltonian system in which u and A have independent 
evolution equations. We begin by defining a phase frequency in the normal velocity frame as 

^+^rn■V<^> = ,^. (4.4) 
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The mass density p and the phase </> are canonically conjugate in the Hamiltonian formulation of the 
Landau two- fluid model. Therefore, = — 6h/5p for a given Hamiltonian h and u = V0 satisfies 

— + v„-Vu+ Vv„f •u = -V— . 4.5 
at op 



The mass density p satisfies the dual equation 



|.V.„v,.).-V.^^, (4.) 

In the last two equations, we see the expected two-cocycle terms for Landau's perfect superfluid, in 
which u: = 0. Here u and p in the superfluid play the roles of 7 and a for the PCFs. The curvature 
u: is nonvanishing now, because of the vortices (defects) that are represented by A. 

Perhaps not surprisingly, the rotational and potential components of the superfluid velocity 
must satisfy similar equations, but the rotational component must be advected by another velocity 
— the vortex line velocity - instead of the normal velocity v„ that advects u. Absorbing all 
gradients into u yields the form of the equation we should expect for A, 

dA 

— + vfxa; = 0. (4.7) 



Taking the difference of the equations for u and A then recovers equation ( [4. 3D as 

v^xa; = — V(v„-u-|- — I with = u — A , (4. 



dt ^ V " ' 5pJ 

in which regularity of the phase (j) allows one to set curl u = 0. It remains to determine vg from the 
Euler-Poincare formulation. Including the additional degree of freedom A allows the vortex lines 
to move relative to both the normal and super components of the fluid, and thereby introduces 
additional reactive forces associated with the momentum of the vortex lines. 

For superfluid ^He with vortices, the momenta conjugate to the velocities v„ and shall be 
our basic dynamical variables. To develop the Euler-Poincare formulation of this problem, we must 
consider a Lagrangian that first of all is invariant under the order parameter group O = U{1). The 
Lagrangian must also be invariant under two types of diffeomorphisms: one corresponding to 
the material motion of the normal fluid G„ and another corresponding to the motion of the vortices 
Gi. Thus we consider a Lagrangian that allows the following direct product of group reductions, 



TGnXV* ><{TO/0))/Gn) X [{TG,xV*)/G, 

^ (gn X {V* X o) g-\t)) X (g, X V*gf{t)) . 

We denote the corresponding dependence in this Lagrangian as 

/ (vn,^, I/, u; v^,n) . (4.9) 

Here and g^ denote the Lie algebras of vector fields associated to the velocities v„ and v^, 
respectively. Also o denotes the Lie algebra of the Abelian gauge group U{1); so o contains v 
and u. The V* in each factor denotes the corresponding advected densities: entropy S advected 
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by v„ and vortex inertial mass n advected by v^. We denote the inverse right actions of the two 
diffeomorphisms as g^^it) and gj^{t). At a given time t, these actions separately map spatial 
variables back to coordinates moving with the normal material and with the vortices, respectively. 

According to the Euler-Poincare action principle, the following dynamical equations are gener- 
ated by this Lagrangian, cf. 



TT- T — = - ad* 



dt 6Vn 
d_6l_ 

d_5l_ 
dt Svi 



and (lOl ) 

61 

div I V. 



51 ^ 51 51 
+ — oS + — ou + ^ 



O Ur. 



5l_ 



+ 



5l_ 



ad* — + — o n 

5v£ 5n 



(4.10) 



These are the Euler-Poincare equations for a superfluid with vortices. The Eulerian kinematic 
equations are, cf. ( 2.61 - 2^ ) 

'dt 
du 

'dt 

dn 

'dt 



-div(5v„) , 
Vn X curl u — V ( U • V 
— div(n v^) . 



If V£ and n are absent, these equations reduce to the equations for a PCF with broken U{1) 
symmetry. The momentum density conjugate to the frequency v is the total mass density given 
by p = —5l/5v, which satisfies the equation above. So far, this is interesting, but standard in 
the present context. However, now a new feature develops because of the physical description of 
superfluids. Physically, nothing is known on the Lagrangian side about the relation of the gauge 
frequency v to the other variables. However, on the Hamiltonian side we know from the Legendre 
transformation that u = —5h/5p. Moreover, the thermodynamic energy on the Hamiltonian side is 
a known function of p. This means we should leave the Lagrangian side to finish determining the 
dynamics for superfluid "^He with vortices. The following Hamiltonian description for this dynamics 
is derived in Holm [2000]. 

Proposition: The dynamics for superfluid ^He with vortices follows from a Lie-Poisson bracket 
whose Hamiltonian matrix separates into two pieces given by 



(4.11) 





' Mi ' 




' Mjdi + djMi 


Sd, 


pdi 


djUi - Uj^i 




" 5h/5Mj ' 
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S 

















5h/5S 


dt 


P 




djp 












5h/5p 




Ui 




Ujdi + Uij 





di 







5h/5uj 



and 



d_ 
di 







' Njdi + djNi no 


n 




djU 



5h/5Nj 
5h/5n 



(4.12) 



^The first of tfiese is tfie ISamiltonian matrix for a PCF witfi broken (7(1) symmetry. Tlie otlier Hamiltonian matrix 
gives the standard semidirect-product Lie-Poisson bracket, without two-cocycles. This combination of Hamihonian 
matrices was first introduced in Holm and Kupershmidt [1987] for superfluid plasmas. 
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where M = 6l/6vn, N = 6l/5vi, p = —5l/5v and the Hamiltonian h is the Legendre transform of 
the Lagrangian I in ^.!\ ) with respect to Vn, V£ and u. 

Corollary: // the Hamiltonian has no explicit spatial dependence, then total momentum conser- 
vation holds as, 

I (M, + AT,) = {M, +N,M = -^, T/ . 
Suppose the Hamiltonian density has dependence h(M., p, S,n,Vs,u, A), in which 

Vg = u — A , A = — N/n and cj = curl . 
Then, the stress tensor T-' is expressed in terms of derivatives of the Hamiltonian as 

rT.k dh f dh / , dh\ \ dh , dh 

= ^i^ir + 7^ + curl — - Vsi,jeraik ^ + 5 P - Aj — - 

9Mfc \dvsk ^ duJkJ '■' duJm dAk 

where the pressure P is defined as 



„ , , dh dh ^dh dh 

Thus, the dependence of the Hamiltonian h on the vorticity u introduces reactive stresses due to 
the motion of the vortices. 

Details of the choice of Hamiltonian, as well as the derivation and interpretation of the explicit 
equations are given in Holm [2000]. Other results include the transformation of these equations 
to a rotating reference frame and the resulting Taylor-Proudman theorem for superfluid ^He with 
vortices. The generalization of this idea to complex fluids with nonabelian broken symmetries will 
be discussed elsewhere. 
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Appendix. External torques and partial Lagrangian reduction. 

The anisotropic dielectric and diamagnetic effects on the director angular momentum due to ex- 
ternal electric and magnetic fields can be restored by adding the torques from equation ( ^.38 ) to 



the right hand side of the second equation in ( p. 52 ). Knowing they can be restored this way, one 



could simply ignore the external torques. However, the major applications of liquid crystals involve 
these torques and their restoration also provides an example of partial Lagrangian reduction. 
This example also briefly recapitulates the procedure of Lagrangian reduction by stages used in the 
remainder of the paper. 
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Partially reduced Lagrange-Poincare equations. Restoring the effects of external torques 
requires that we consider Hamilton's principle for a Lagrangian that still retains its dependence on 
the director n, 



S' = j dt j d^X £(x, J,n,i/,Vn) 



(4.13) 



In this case, the partial reduction of the Lagrangian dependence from h to i/ = n x h proceeds as 
follows. The variation of = n x n gives 



(Ji/ = (n X 5n)' — 2n X 5n = (n X Sn)' — 2i/ x (n x 5n) . 
Moreover, we have the relation 

(^n • A = (n X 5i\) • (n x A) = |np(5n • A — (n • (^n)(n • A) , 



(4.14) 



(4.15) 



for any vector A, since |np = 1, which implies that n • (5n = 0. Substituting the identity ( 4.14]) for 
5u and the relation ([4.15|) into Hamilton's principle implies the following replacements, 



\ an 

[it 

Van 



(n X 5n) 



d, 



(n X 5ri) • ( n X ^) 
an/ 



(4.16) 
(4.17) 



Varying the action S' in the fields x, n and v at fixed material position X and time t now gives 
5S' = - I dt I SXI 5xr 



dx 

+ (n X 5n) 
— (n X Sn) 



d£\ d dC d f . dC 

+ Jt^— - J- I J n, 



dxp dJ dxm V 9n.i 



_2u X — 



nx(—- jJL(j-^^ 
V 5n dxm V 9n,m 



(4.18) 



with the same natural (homogeneous) boundary conditions as before. Consequently, the action 
principle 5S' = yields the following partially reduced Lagrange-Poincare equations for 
liquid crystals. 



5xp : 



dC\ d dC d ( . dC 

+ Jt^— - J- — I J' n, 



dx„ J dxp dJ dxn 
dC\ . dC 



dn , 



0, 



(4.19) 



nxJn: \ — \ - 2u x nx - J 

ov ov \ on ox 



dC 



d 



J 



-1 



dC 

dnm. 



0. 



One may compare these with the more completely reduced Lagrange-Poincare equations for liquid 
crystals in equations ( 2.52| ). The explicit dependence of the Lagrangian on the director n introduces 
torques and stresses not seen for Lagrangians depending only on i/ and 7 = n x tin. 
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Partially reduced Euler-Poincare equations. We are dealing with Hamilton's principle for 
a Lagrangian in the class 

S' = j dt j d^x i{u,D,n,v,Vn) , (4.20) 

in terms of the Lagrangian density £ given by 

i{u, D, n, jy, Vn) d^x = (4.21) 
c(xg-\t),Jg-\t),ng-\t),,yg-\t),Vng-\t))(d^Xg-\t)). 



(4.22) 
(4.23) 
(4.24) 
(4.25) 



The variations of the Eulerian fluid quantities are computed from their definitions to be, 



6uj = 
5D 

5u 

5n 



dt oxk 
dDrjj 
dxj ' 

as ^ ^ as 

dt dxm 



dxk ' 



2u X — rjr, 



du 

dXm. 



dn 

dxk 



where Ti{x,t) = (n x 5n){X,i)g'^ {t) and = 5gg~^{t). 

We compute the variation of the action ( |4.20| ) in Eulerian variables at fixed time t and spatial 
position X as, cf. equation ( 2.65| ), 



5S 



dtld^x 

3 



di , 6£ 6£ ^ 6£ ' 

^Suj + Y7^6D + — -5v + — -6n 
ouj oD ov on 



dtid xirjj 



d 5£ _ _5^duk _ _d_ 
dt 6ui 



6uk dxj dxk \Suj 



+ D 



dxj 6D 



5£ du 
5iy dxj 



6£ dn 
6n dxj J 



(4.26) 



+ S 



d_6£_ 
dtJu 

d 



d 



dXm ^SU 



-Ur, 



6£ 



+ 2v X — + nx 

ou on 



+ 



dt 



5u. 



+ 



d 



dXr. 



+ S 



5£ 



6£_ 
5^ 



6u 



where we have substituted the variational expressions ( [4.22D -( p35 ) and integrated by parts. Hence, 
we obtain the partially reduced Euler-Poincare equations for liquid crystals, 



d b£ 
dt 6ui 



6£ duk 



d / 6£ 



S- - — 
■ dt6u 



5uk dxj dxk ^ Suj 
6£ du 
6u dxj 

d 



Uk) + D 



d S£ 
dxj 6D 



6£_ dn_ 
<5n dxj ' 



5£ \ 6£ 5£ 

a V -^^rn + 2U X —+ nX —. 

dXm \du / du on 



(4.27) 



(4.28) 
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One may compare these with the more completely reduced Euler-Poincare equations for liquid 
crystals in equations (|2.66| ) and ( 2.67| ). Again ones sees the torques and stresses generated by the 
explicit dependence of the Lagrangian on the director n. 



Partially reduced Hamiltonian dynamics of liquid crystals The Euler-Lagrange-Poincare 
formulation of liquid crystal dynamics obtained so far allows passage to the corresponding Hamilto- 
nian formulation via the following Legendre transformation of the partially reduced Lagrangian 
i in the velocities u and v, in the Eulerian fluid description, 



mi = — , a = — , h{m, D, a, n) 
OUi ou 

Accordingly, one computes the derivatives of h as 

5h 



rriiUi + a • v — i{u, D, v, n). 



6h 
6^ 



V . 



6h_ 

Id 



Sh 
6^ 



(4.29) 



(4.30) 



5mi 6a 6D 6D 

Consequently, the Euler-Poincare equations ( p. 66 ) - ( 4.2^ ) and the auxiliary kinematic equations 



(2.61) - (2.62) for liquid crystal dynamics in the Eulerian description imply the following equations, 
for the Legendre-transformed variables, {m,D,a,n), 



drrii 



dD 
'dt 
dn 
'dt 
da^ 
'dt 



d 6h 



dxi 5m j 



_d_ 

dxi 



mi 



/ dn 

\dxi 



D 



d 

dxj 
dn 



dxj J 6m j 



6h 
6n 

5h 
6mj 
6h 



6h 
5m j 
d 6h 
dxj 5a 



dxj 5D 



n X 



6h 
5^ 



_d_ 

dxj 



5h 
6m i 



n X 



5h 
5^ 



2a X 



5h 
5^' 



(4.31) 
(4.32) 
(4.33) 
(4.34) 



These equations are Hamiltonian. Assembling the liquid crystal equations ( [4.31| ) - ( 4.34 ) into the 
Hamiltonian form ( p.79| ) gives, 





mi 




d 


D 




di 


n 






a 





m 



jdi + dj 
d,D 

dja 



mi 



Ddi 







-Hi 




nx 



a • d^ 


nx 
2a X 



5h/5mj 
5h/5D 
5h/5n 
5h/5a 



(4.35) 



One may compare this with the more completely reduced Hamiltonian matrix form for liquid crystals 
in equations ( p. 81 ). The Jacobi identity for the Poisson bracket defined by this Hamiltonian matrix 
is guaranteed by associating it to the dual of a semidirect-product Lie algebra that, in this case, 
has no two-cocycles. The two-cocycles are generated by the further transformation to 7 = n x dn. 
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